Postgraduate Lectures 

1. Time Independent Green's functions 
a) Formalism 

Green's functions are denned as the solutions of the inhomogeneous equation: 

[z-C(r)]G(ry-z) = Ur-r') (1.1) 

where z = X + is is a complex variable and £(r) is a time-independent, linear, Hermitian operator with eigenfunctions 
<j)„(r) and 

C{r)Ur) = KUr) (1.2) 

(])„(r) and G{r,r';z) satisfy the same boundary conditions on the surface of the domain Q.. The set {(])„} can be 
considered orthonormal and complete i.e. 

/ <t£(r)<|> m (r) dr = 8 nm (1.3a) 

^„(r)C(r') = o(r-r') (1.3ft) 

n 

Note that the sum J2 and integral / operations are more or less interchangeable depending on the nature of the operator 
C, continuum or tight binding, and its spectrum, continuous or discrete. It is useful to consider an abstract vector space 
and use Dirac's notation 

4>»M s <rM>»> 

b{r-r')C{r) = {r\L\r') (1.4/?) 

G{ry-z) = {r\Q{z)\r') (1.4c) 

{r\r') = 8(r-r) (I Ad) 
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dr\r)(r\ = l (lAe) 

(z-£)G(z) = l (1-4/) 
£|^) = X„|^) (lAg) 

<<M<|>m> = 8nm ( 1 -4/z) 

X>«><4>«| = 1 (1-40 

We can revert to the /--representation by taking the (r\, |r'> matrix element, thus 

<r|fe-£)e(z)|r') = <r|l|r') = <r|r') = 8(r-r') (1.5a) 
zG(r, r'\ z) - (r\ CQ{z) \r) = zG(r, r'; z) - J dr" (r\ £ \r") {r" \ g(z) |r'> 

= zG(r, r';z)- C{r)G{r, r';z) = 8(r - r') (1 .5b) 

If z + {X„} then we can write 

G(z)=^- T (1.6) 

z- £ 

Multiplying (1.6) by (1.4i) and using F(£) |(j)„) = F(k n ) |(j)„) gives 

G(z) = ^J2\^^»\ 

n 

= X] ~T£ I *!>«> <*t>« I 
_ ^ | *(>«> <<t>« | 



(1.7) 



G(ry;z) = J2 



z-K 

WWn(r) 

z-K 
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Since £ is hermitian, all its eigenvalues {X n } are real. Hence Q(z) is an analytic function everywhere on the 
complex plain except at those points which correspond to the eigenvalues of £. Q(z) has poles at the the discrete 
eigenvalues of £. The converse is also true. 

If the spectrum is continuous, i.e. there is a range of values of X such that every value in that range is an eigenvalue 
of £, Q(X) is not well defined. If we try to define it by a limiting procedure we find that the side limits lim.,_>o + G(X±is) 
exist but are different from one another. This is a branch cut of Q(z) along the real z-axis. 

In disordered systems there are other possibilities besides discrete and continuous spectra (see later). 

We can define 2 Green's functions as follows 

G + (k)= Mm Q(X + is) (1.8a) 

Q~(X)= Mm Q(X- is) (1.8b) 

From (1.7) we obtain the results 

G*(r,r';z) = G(r',r;z*) (1.9a) 

G-(r,r'; X) = [G + (r',r; X)] * (1.%) 

ReG-(r,r;X) = ReG + (r,r;X) (1.9c) 

lmG~(r,r,X) = -lmG + (r,r;X) (l.9d) 

lim — !— = P - =f ix8(x) (1.10) 

v^o+ x+iy x 

g(X) = g + (X)-g-(X) (l.i i«) 

G(r, r';X) = -2%i £ b(X - X n )WW n (r') (1.1 lb) 

n 

Integrating the diagonal elements over r gives 

Gfy.r; X) = P £ T in £ 5(X - 5l„)W)4>» d-12fl) 

n n n 

TvG ± (X) = Pj2x^Y^ i% J2 5{X - X " ) {Ll2b) 

n " n 

J2„ 8(A, - X n ) is the density of states (DOS) at X, N(X); N(X) dX is the number of states in the interval [X, X + dX]. The 
quantity 

p(r;k)^5a-^„(r)4)>) (1.13) 

n 

is the DOS per unit volume, and N(X) = J p(r; X) dr Hence 

p(r;X)=T-lmG ± (r,r;X) 
% 

= ^-.G(r,r;X) (1.14a) 
2%i 

N(X) = T - Im Tr g\X) (l.Ub) 



Using the identity 



we can define 



K 



Another important relationship is 



G(r,r';z) = J2 



§ n (r)§* n (r') 
z-X„ 



-L 



dX^-K)^^ (LIS) 

n 

G(r,r';X) 



i 



z-X 



dX- 



2n z-X 
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In particular 

G(r,r;z) = jf"dX^f£ d-16) 

Thus the real and imaginary parts of £ are not independent of one another. (1.16) is an example of a Kramers-Kronig 
relation. One more important property is 

dG(r,r ,z) = _ J ^„ G{r y u ) G {r"y u ) (Ilia) 
WWrSr") i(r")C(r') 



« m 



E 



(z-Kf 



cfe Z — ' z - X„ 

n 

= -{r\(X-jr,y 2 \r) 



dG(r,r;X) , , 



dA, 

<0 X±{X n } (1.17b) 
knowledge of the Green's function permits us to solve the general inhomogeneous equation 

[z -£(#•)] M (r)=/(r) (1.18) 
where u(r) satisfies the same boundary conditions as G and (j). By operating with Q on (1.18) we obtain the solution 

u(r) = J G(r,r';z)f(r')dr' z*{U (1.19a) 

GHr.r'iDfir^dr' + ^r) z = {X n } (1-1%) 



where X n belongs to the branch cut, the continuous spectrum, of £, and (j)(r) is a solution of the corresponding 
homogeneous equation. If z coincides with a discrete eigenvalue of £ , say X n , there is no solution of ( 1 . 1 8) unless /(r) 
is orthogonal to all eigenfunctions (j)(r) corresponding to X n . 

If u(r) describes the response of a system to a source f(r), then G(r, r) describes the response of the same system 
to a unit point source located at r' . Note that because of ( 1 .9) this response is essentially the same as that at r' due to a 
point source at r. (19) implies that the response to f{r) can be expressed as the sum of the responses to point sources 
distributed according to f(r). 
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b) Example 1 

Consider the case £(r) = -V 2 and the domain Q. extends over the whole 3 dimensional space. Then the eigenfunctions 
of £ are 



(r\k) = -^=e ikr (1.20a) 



with eigenvalues 

X n = k z (1.20b) 

where k is real. The spectrum is continuous, from to +°°. 

The Green's function can be obtained either by solving the defining equation 

(z+V 2 )G(r,r';z) = 8(r-r') (1.21) 

or from the solution of the homogeneous eigenvalue problem using (1.7) 

<#•!*> <*|/> 



* z ~ k2 



(2%) d z-k 2 



1 



(1.22) 



where d is the dimensionality. 

Consider the 3-dimensional case d=3. If p = r - r' and 9 is the angle between k and p we can write ( 1 .22) as 



1 t°° k 2 Mr f n 

v-^w-L t^-L **** 



cos 6 

11 v/c 

k 2 dk e +ikp - e~ ikp 



z-k 2 ikp 
1 r°° ke ikp 

dk 



(1.23) 



4in 2 p J_„ z-k 2 

the integration path can be closed in the upper half plane. Unless z is real and nonnegative, one of the poles, +y/z, has 
a positive imaginary part and the other, —y/z, has a negative imaginary part and lies outside the integration contour. 
Hence 

G(r,r';z)= , jr ImVz>0 (1.24) 

4jc|r-r| 

For z = X and X > G is not well defined, but the side limits are 

±«V%-/| 

G ± (r,r';X) = , -r VX>0, X>0 (1.25) 

4ti r - r 



For z = X and X < we obtain 



e -Vty-r>\ 



G(r,r';X) = , -r VX>0, X>0 (1.26) 

4n\r-r\ 



and for z = we have 

G(r,r';0)= — 1 (1.27) 
4jc r — r 

which is the Green's function corresponding to Laplace's equation 

V 2 G(r,r ;0) = 8(r-r) (1.28) 
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By applying (1.19b) we can find the general solution of Poisson's equation 

V 2 V(r) = -— p(r) (1.29a) 
£o 



V(r) = / G(r, r'; z) — p(r') dr' + C 
J £o 

1 / p(r') dr' 



£o J \r-r 



+ C (1.29/7) 



The constant has been added since the general eigenfunction of-V 2 with eigenvalue is a constant. (1.29) is the basic 
result of electrostatics. 

c) Example 2 

Now consider the 1 -dimensional case 

z + ^] G(x,x';z) = 8(x-x) (1.30) 



dx 2 . 

For x + x the solution must be that of the homogeneous equation 

G = Ae +i ^~ x + Be- i ^ x (1.31) 

With Im i/i > only the + solution is finite as x - x' —> +°° and only the - solution for x - x' — > By integrating 
(1.30) we find 

G(x - ,x';z) = G(x + ,x';z) (1 .32a) 

£).-£) 

/ x-x- / x=x + 

We can thus determine the coefficients A and B to obtain 

giVi|*-x'| 

G(x,x';z)= - n Imv / i>0 (1.33) 
2itJz 



For z = A. > we have for the side limits 



G ± (x,x;X) = T-^i=e ±iVx \ x - x '\ X>0 VX>0 (1.34) 
2vA- 



whereas for z = - I X\ we obtain 



G(x,x';-\X\) = l — e -^M\ x - x '\ X<0 J\X\>0 (1.35) 

2^/\X\ 



and for the 1-D Laplace's equation, z = 0, 

G(x,x';0) = i Ix-x'j + const (1.36) 
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2. Time-Dependent Green's Functions 
a) First-Order Case 

The Green's function g(r,r , t - f) associated with the 1st order (in time) P.D.E's 



is defined by 



cdt 
i_d_ 

cdt 



i d 
c at 



-C(r) 



W. = 
Mf(r,t)=f(r,t) 



g(r,r , ,t,f) = 8(r-r')8(t-f) 



(2.1a) 
(2.1b) 

(2.2) 



subject to the same boundary conditions as (j) and \\t. For the moment assume that c is a positive constant. The operator 
C(r) is as in section 1. 

Expressing g(r, r , x), where x = t - t' via its Fourier transform we obtain 



g«i» = Hr-r) 

g m= G n 



(23a) 
(2.3b) 
(2.3c) 



where G(z) is the Green's function discussed in section 1. As before it may be necessary to use a limiting procedure 
to define g(z). 



In analogy with the previous procedure we can also define 

g(z) = g + (x) - g~(x) 



(2.4) 



(2.5) 



Note g is integrated along a line just above the real axis and back just below. This contour thus contains all the 
eigenvalues of C(r). However g ± (l) requires its integral to be closed around either the upper or the lower half plane. 
This is chosen such that the exponential becomes negligible well away from the real axis. Thus the appropriate choice 
is the lower (upper) half plane for x > (x < 0). 
The g's obey the relationships 



^=±0(±x)|(x) 
g-(rS,x)= [g + r',r-z)\ 



(2.6a) 
(2.6b) 



We can write g(z) as 



dco' 
~2% ( 



= -2ni jT 1^ £ 5 (7 " *») UrW n (r) 
= -icY j e- icK 'Um* n (r') 

n 

The corresponding operator g(x) is then 



(2.7) 



g(x) = -ic^V 



icX n i 



0«> <<|>J 



(2.8) 



= -ice 



-iCx 
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where 

U{t - f) = e - icC{, - f) = -g(f - f) (2.9) 
c 

is the time-evolution operator or a propagator because 

\§(t)) = U(t-t')\W)) (2.10) 

Note also that 

g(r,r' ,0) =-ic8(r-r) (2.11a) 

u(h - r 2 ) = m x - hmh - 1 2 ) (2.iib) 

By combining (2.9) and (2.10) we obtain 

4>(r, t) = l - J g(r, r', t - t')<?(r\ f ) dr' (2.12) 
The solution of the inhomogeneous equation, \|/(r, f), can be expressed as 

\|/(r , f) = (j)(r, t) + J dr' df g + {r, r , t - t')f{r', t') (2.13a) 

= <t>(r, t) + y dr' y At' g(r, r\ t - t')f(r', f) (2. Ub) 

This is best proved by substitution of (2.13) in (2.1) and using (2.2). If we had used g~{l) instead of g + {i) the resulting 
\|/(r, t) would have satisfied (2.1b). We have excluded this solution using the physical argument that the response of 
the system at a time t depends only on what the source was in the past ( f < t). 

b) Example 

Consider again the case C = -V 2 . It suffices to calculate g(%). Using (2.7) and (j)„(r) = exp(;'£ • r)/y/Q. and X n = k 2 we 
get 

. . p 'kir-r') 

g { r,r'^)=-icY J —^ r ^ k ' 

k (2.14) 

= - ic / _^J*-P-« 2 T 

J (2%) d 

where p=r-r' and d is the dimensionality. Using k ■ p = J2t=i kPi an< ^ ^ = we can rewrite (2. 14) as 



g(r,r',x) = -icf[f "§e* pH * 2x 

i=l ~°° 

ni£d 1 Fk~ 
?TT \/ ;7-T 



= -ic | | / —e^e- 1 '^ 2 

(2.15) 



271 V 'CT 

(=1 



1 



4tc/ct 



d 



•JL. 

e 4 " 



where ^/nlicx is the square root with a positive real part. (2.15) and (2.12) allow us to study the time evolution of a 
wave-packet. 

If we choose c = —iK where K is a positive constant, and £ = V 2 we obtain the diffusion equation -3(j)/K3f + V 2 <j) = 0. 
Again for causality reasons we need g + (x) only. It is more convenient to consider the Laplace transform, rather than 
the Fourier transform. The final result, obtained by substituting a— > - z'K in (2.15), is 

1 \* f-(r-r'Y 



Note that g + >— > - K8(r - r') as t - ft— >0 + . As t increases, g + describes the diffusion of a local initial pulse. Note that the 
average diffusion range increases as \/t-t' '. 
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c) Second-Order Case 

The second-order differential equation gives 



■Or) 



' c 2 dt 2 



c 2 dt 2 



dt 2 



W> = 

\lf{r,t)=f(r,t) 



g(r,r\t-t') = b{r-r')b{t-t') 



(2.17a) 
(2.17/j) 
(2.17c) 



subject to the same boundary conditions on (j), \\f, and g. As before we express g(x) as a Fourier transform of g(co) 

' +0 ° dco' 



V«') 



S(co) = G 



co- 



(2.18a) 
(2.18fe) 



Since G(z) has all its singularities on the real z-axis, g(co) is analytic everywhere on the complex CO plane except 
on the real or imaginary axes. However on physical grounds we can restrict ourselves to the case where all the 
eigenfrequencies co„ are real. Then the singularities of G(z) are all on the positive real z-semiaxis. 

Again near a branch cut it may be necessary to use a limiting procedure to define g(0)). There are 4 cases 



g + (co') = lim G —r- + is 



.5^0+ 



co- 



co 



lim g (co' + we(co')) 



.5^0+ 



/(CO') = lim G 

.s->0+ 



g R (a>') = lim G ( —z — h /se(co') ] = lim g (co' + is) 



g (co') = lim G —j- - is 



CO 



co- 



ise( co') = lim g (co' - is) 



= G~ 



co - 



= lim g (co' - ise( co')) 

,!->0+ 



where 



L— 1 if; 



(2.19a) 
(2.1%) 
(2.19c) 

(2.19a 1 ) 
(2.20) 



ifx>0 
? x<0 

and co' is real. These are referred to as the causal (a), retarded (b), and advanced (c) Green's functions. The fourth 
case is not independent and can be derived from g R + g A = g + + g~. 
We can also define 3 different g" s as 

A 



8 = 8 -I 
t = 8 + -8 R 



8 = 8 R -8 A = 8 > ~8 < 



(2.21a) 
(2.21/7) 
(2.21c) 



Note that the g's satisfy the homogeneous equation (2.17a) rather than (2.17c). 

As before by noting that the path for the various g's can be closed in the upper (or lower) half plane when x is 
larger (smaller) than zero we have 



g + (t) = 6(T)f(T) + 0(-X)r(T) 

g R (x) = 0(x)f (x) 
/(x) = -0(-x)|(x) 

f(T) = -6(T)r(T)-0(-T)f(T) 



(2.22a) 
(2.22/7) 
(2.22c) 
(2.22d) 
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It can also be shown 



Jo 

" jf ^ e ^ 2TO ' £ W^>'> 6 " X «) (2-23) 



2 V 



where ^ 0. Similarly 



and also 



* 2 4- vx; 



t(r,r',x)=-[g > (r',r,x)]* (2.25a) 

WW) „ ; „ ^ ^ (2 25fo) 



g(r, r , x) = -c j= — sin [c y/X n zj 



sin 



g(x) = -c -j= in operator form (2.25c) 



Consider now the expression 



|4>(f)> = [m-t') j(f.(fO) + g(t-t') (2.26a) 
4>(r, t)=~jdr' g(r, r', t - f)^r\ f) 

- ^ J dr' g[r, r\ t - t')4?{r\ f) (2.26b) 

where the dot denotes differentiation with respect to t and (j>(f') is d<S?l&t for t = f. Since g(t-f) satisfies the homogeneous 
equation (2.17a) so does (j)(r, f). Also (j)(f) — > §(f) and (j>(f) — > (j)(f') as t — > f. Thus this determines the solution of 
(2.17a) for time t in terms of §(f) and (j>(f') at time t'. Then 

\|/(r, t) = (j)(r, f) + y dr' At' g R (r, r, t - t')f(r\ f) (2.27a) 
= (j)(r, t) + J J g(r, r, t - t')f(r', f) (2.21b) 

satisfies the inhomogeneous equation (2.17b) where (j)(r, t) is the general solution of (2.17a). g R is used because the 
response can only depend on the source at times t' < t. 



9 



Transport & Disorder 



Chapter 3. 



c) Example (2nd Order) 

Consider again the case £ = -V 2 in 3-D for which G ± (X) are given by (1.25). Substituting in (2.23) we get 

g>(p, x) = -— f ^. e ~ im ( e +i< °P /e - e- i( °P /c ) (2.28a) 
4%p Jo 2% v > 

|<(p, x) = I ^V ,<0T ( e +,a, P /e - e- i( °P /c ) (2.28b) 

where p = r - r . Subtracting we get 

ff(p, x) = — / — (e +koplc - e- iaplc ) (e +i(m - e- im ) 
4np Jo 2ti v y v ' 



_ /" dco / !t0( p/ c+T ) _ e ,(o( P / e -T)^ 

4;rp ( c + X ) ^ ( c T ) 
= 4^p~ [ 5 (P + CT )- § (P- CT )] 



Using (2.29) and (2.22) we obtain 



4jc |r - r 

Then the solution of the inhomogeneous wave equation 



is from (2.27) 



V|/(r, f) = <|)(r, f) - — / dr' df' 8 ( £ - (t - f)) 

4% J V c / p 



which is the basic result of electromagnetic theory. 



(2.29) 



g s (r,r', f - = i ttS (|r - r'l - c(t - t')) (2.30a) 

4jc r — r 

sV.r', f - = p — 7T§ (|r - r' | + c(t - t')) (2.30a) 



(2.31) 
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3. Application to Free Quantum-Mechanical Particles 
a) General Relations 

The one-particle time-independent Schrodinger equation and its Green's function have the form 

[E-7i(r)]W) = (3.1a) 
[E-Ti(r)]G(r,r';E)=8(r-r') (3.1b) 

where Ti(r) is the Hamiltonian in the /--representation and G(r,r; E) satisfies the same boundary conditions as \|/(r). 
From section 1 . it is clear that 

C(r) -> Ti(r) {3.2a) 

X^E (3.2b) 

X+is = z^> z = E+is (3.2c) 

K -> E n (3. 2d) 

4>»(r) -> Mr) (3.2e) 

Thus G is related to the eigenvalues, E m and eigenfunctions, (j)„, of Ti by 

C(ry: z) = ^J0p. (3.3«) 

The singularities of Q(z) are on the real axis and 

1) The poles of Q(z) coincide with the discrete eigenvalues of Ti, and vice versa. 

2) The residue at each pole E n of G(r,r';z) is J2i^i( r )^K r ) where J2i i s over tne fn degenerate eigenstates corre- 
sponding to E„. 

3) The degeneracy /„ can be found by integrating the residue of the diagonal matrix element G(r, r;E„) over r i.e., 



f n = J d 3 r Res G(r,r;E n ) 



(3.4) 
= TrRes6(£„) 

For a non-degenerate state /„ = 1 and 

UrW„(r') = R^G(r,r';E n ) (3.5a) 
\Ur) | = >/[ResG(r, #■;£„)] (3.5ft) 

M ., f ResG(r,0;£„) | 

PnW = In < > (3.5c) 

I VfRes G(r, r ; E n ) ■ Res G(0, 0; £„)] J 

where /?„(r) is the phase of (j)(r) (assuming p„(0) = 0). 

4) The branch cuts of Q(z) along the real z-axis coincide with the continuous spectrum of Ti and vice versa. 

5) The density of states per unit volume is given by 

p(r;E) = T-lmG ± (r,r;E) (3.6) 

K 



6) The density of states N(E) is given by 

N(E)= I drp(r;E) 



= / drp(r ;j 

J (3.7) 
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7) The time-dependent Schrodinger equation gives 

.a 



\y(t)) = U(t-to) |\|/(f )> 

U(t-t ) = exp I -i J 

= ifigit - t Q ) 

= iTl [ + °° ^ e -i«xt-to)Q (m 

\|/(r, f) = ih J g(r, r, t - t )y(r', t ) dr' 



b) The Free-Particle Case, % a = p 2 /2m. 

We denote the free-particle Hamiltonian by Tio where 



° 2m 2m 



and Go(z) is the corresponding Green's function 



z+^V?) G„(r.r':,-.) = 5(r-r') 
Hence 

where G{r,r; z) is the Green's function corresponding to the operator £ = -V 2 . Thus we have for Go 
3-D Case 

m <T*°l r - r 1 

G Q (r,r';E) = T1 w E<0 

2th 1 \r-r\ 

m e ±ik "\ r ~ r '\ 

Gt(r,r';E) = T -r- E>0 

27tfi 2 r-r' 



<7o(z) has a branch cut along the +ve £-axis corresponding to the continuous spectrum of Tio. Then 

G(r, r';E) = G + (r, r';E)- G~(r, r'\ E) 

m sin (kn I r-r' \) 
= -2ni - y l ]) ®(E) 

2nti \ r - r \ 

G Q (r,r';E) 



po(r;£) = 



-271/ 

2n 2 h 2 

V2tc 2 7z 3 



1 Ti = h/2n, where h is Planck's constant 
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1-D Case 



G Q (x,x; E) = -T-e^M E < (3.14a) 
Ti 2 k () 

GUx,x; E) = T -^- e ±( *°M E<0 (3.14b) 

h k a 



Gq(x,x';E) = -2ni&(E) — = — cos (ko Ix-x'l) (3.14c) 



p Q (x; E) 



m 
' Kh 2 k 
Gq(x, x; E) 



0(F) 



2ni 
m 



nh 2 k Q 



Note that the behaviour of p(r; E), and hence of G(E) and G ± (E), at the band edge E — > + depends crucially on the 
dimensionality. In general 

lim 0(E) ~ (E-Eb)™ (3.15a) 

{ const, in 3-D 

-ln(z) in2-D (3A5b) 

4. Perturbation Theory 
a) Time-Independent Case 

Suppose that the one-particle Hamiltonian, Ti, can be separated into an unperturbed part, Tio, and a perturbation, Ti i 
where 

n = n Q +ni (4.i) 



Assuming we can solve for the Green's function Go of Tio, how can we find Q of Til Noting that 

G (z) = (z-Ti Q r l (4.2a) 
G(z) = (z-TiT 1 (4.2b) 



we proceed as follows 

Q = [z-Ti a -Ti l Y l 

W-Ti,]Q = l ( 4 - 3 ) 

[I-GoHi\G = Go 

G = G0 + G0H1G 

This is often called Dyson 's equation. It is the starting point for almost everything that follows. Note that (4.3) can be 
expanded to give 

G = Go + GoTi l Go + GoTi 1G0 + . . . (4.4a) 

= Go + (Go + GoTiiGo+-..)'H l Go (4.4b) 

= Go + GTi 1 Go (4.4c) 

In the /--representation (4.3) becomes 

G(r,r';z) = G (r,r';z) + J dr\ dr 2 G a (r,r x ; z)Hi(ri,r 2 )G(r 2 ,r'; z) (4.5) 

13 
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Often H\(r\,r 2 ) is a local potential which can be written 8(ri —r^Vifi); then (4.5) becomes 

G(r,r';z) = G Q (r,r';z) + J dr, Go(r,n ; z)y(r 1 )G(r 1> r'; z) (4.6) 

which is a linear inhomogeneous integral equation with a kernal of the form Go{r,r\ ; z)V(ri). 
In the ^-representation we can write 

G(k, k';z) = G Q (k, k';z) + J2 Go(*. *i ! zW\ *2)G(* 2> z ) (4.7) 



where we have used the relations 



e ft-r 



W = 7H (4 ' 8fl) 



d£ 

, (2J0 5 



(4.8fe) 



Next we introduce the t-matrix T(z) which can be defined via 

S(z) = So(z) + So(z)T(z)So(z) (4.9) 

By comparing (4.9) with (4.3) and (4.4) we see that 

T = ~H\ + H\ (Go + GqH\Go+ ...)H\ 

= %i+'HiG'Hi (4.9a) 
= H\+H\Qq(J-L\ +H\Gq'H \ + ...) 

= %i+%iGoT {A. 9b) 
= %\ +{T~L \ +Ti iGqH i + ■ • •) GaH i 

= 7ii +TG Q n l (4.9c) 

Clearly all the reservations about the behaviour of G (z) close to the real axis also apply to T(z). Thus T(z) is undefined 
where z is equal to a discrete eigenvalue of Ti, and we must define T ± (E) close to a branch cut (a continuum). 
As before (4.9) will become an integral equation in either r or k representations. Thus 

T(k, k';z) = Hi (k, k') + Y J Hdk,k l )G (kuk 2 ;z) T(k 2 , k';z) (4.1 0) 

kik 2 

H x {k,k')={k\ni \ k ')=^ j drdr'e- ikr+ik ' r 'n 1 (r,r) 

G Q (k u k 2 ;z) = <*i|Go(z) l*2> = ^ y ^ dr 2e -*'"-' + * 2 "- 2 G (r 1> r 2 ;z) 

2X*.*';z)s<ik|T(z) |*0= ^ y drdr'e-*"' + * v r(r,r / ;z) 

In the usual case when Ji\{r,r) = 8(r-r')V(r) and Go(r,r';z) is a function of r — r' only. This simplifies to 

Hi(k,k') = V{k-k')^ 

V(q) = J dr V(r)<T !qr q = k - k' 
Go(k 1 ,k 2 ;z) = h lk2 G () (ki;z) (4.11) 
G Q (kuz) = j dpG Q (p;z)e- ik ' p p=r-r' 

T(k,k';z) = V(k-k') + j -^Vik-k^GoikuzlTXkuk'iz) 
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where T'(k,k';z) = £lT(k,k';z). 

From G(z) (or T(z)) we can find the discrete eigenvalues of Ti and the corresponding eigenvectors. In the 
continuous case, only the DOS can be obtained directly from Q{z). However we can write the time-independent 
Schrodinger equation in the form 

(£-Wo)|v) = ^i|¥> (4-12) 
which can be considered as an inhomogeneous equation whose solution is 

|¥ ± ) = \$> + Gt(E)Hi Jv 1 ) (4.13) 

where |(j)) is the general solution of (E - Ho) = and i)/ 1 are associated with Qq respectively. In the /--representation 
this becomes an integral equation 

x V ± (r) = m + j dr 1 dr 2 G^(r,r 1 ;£)// 1 (r 1 ,r 2 )\|/ ± (r 2 ) (4.14a) 
= 4>(r)+ J 'dnGfcr.r^fiMrOxifVd (4.Ub) 

where (4.14b) applies in the usual case of a local potential. (4.14) is the Lippmann-Schwinger equation. If E does not 
belong to the spectrum of Ti (4.14) becomes homogeneous, (j)(r) = 0. 



/ 



y(r)= dnGofr.r^nr,)^,) (4.15) 



By iterating (4.14) and using (4.9a) we obtain 

\v t ) = \*) + Gt\V> + GpiiGpii |<t>>+- 



(4.16) 



= \if) + gp + -\if) 

By multiplying (4.4) from the left or right by Ti \ we have 

Ti l G = rg Q or GTii=G T (4.17) 

which when substituted into (4.16) gives 

\V t ) = \V> + G ± ni\Q) (4.18) 

Comparing (4.16) with (4.13) gives 

T ± U) = W 1 j ¥ ± ) (4.19) 
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b) Time-Dependent Case 

Using similar procedures on the time-dependent Schrodinger equation 



iH di ~ n °) I V(0) = ni(t) ' ¥(f)) (4 ' 20) 
where Hi(t) may be time-dependent. From the general solution of an inhomogeneous equation we have 

/+oo 
&t'gl(j-t')7ii(t')\y)f\t')) (4.21) 

On physical grounds we keep the | \|/ + (f)) solution because it is causal; the effects of Ti i (f) appear after Ti i (f) has been 
applied. 

Iterating (4.21) gives 



|¥ + (0) = |4>M>+ J Ahglit-h)nm \Wx)) 

+ / df 1 df 2 gS(f-fi)Wi(fi)gS(fi-f 2 )Wi(f2)|(l)(f2))+. 



(4.22) 



Assuming Ti\{t) = for t < fo and J(j)o(f)) is an eigenfunction of Ho, say |(j)„), we have 

f E n (t-t )\ , 
|(j)(/)) = exp I -i J |0„> 

=ex P p° ( ;- fo) ) ^ 

= iHo(t-to) |<(>«> 
Under these initial conditions (4.22) can be written as 

\y + (t)) = A(t,t Q )\<b„) (4.24a) 
A(t,t ) = ifii (t-to) + ih / dfi g (r- fi)Wi(fi)g (fi - fo) 

+ m I dtidt2i (t-h)ni(h)g^(h-t 2 )ni(t2)i (t2-to) + ... 

J t 

(4.24b) 

The probability amplitude for a transition from state |(j)„) to the state |(j) m ) as a result of Tii(t) acting during an 
interval [fo, f] can be calculated from (4.24) as 



(E t + E t{) \ 
~i M n " ° j [<<l>ml<|>n> 



x exp I -i— — I |q>„> + ... 
We define the operator S(t, fo) as follows 



dfi <<(> m | exp I i—^- Wi(fi)exp ( -/— ^ ] |(j)„) 



d*i df 2 <<|>J exp I i^p- ] Wi(fi)gJ(fi - f 2 )Wi(f 2 ) 



(4.25) 



5(/./„) exp ( i^j A(t,t )exp (-i^ ] (4 26) 
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so that ((j) m | S(t, to) |<t> n ) is the quantity in parenthesis in (4.25) which can be rewritten as 

(ty m \S{t, f ) |(j)„) = 8 m „ + / d/i e iaWl <(j) m | Wi(f) |(j)„) 



+ - /'df,df 2 / ^^.(<o„-<o V r 2 «o-<o„) (4-27) 



where co„ = £„/^ and (d mn = C0 m - C0„. 

(4.27) implies that the probability amplitude of a transition between two different states as a result of Ti i (f) acting 
during an infinite period, -°° — > +°° is 



<4>m|<S|4>m>=^ / d/,f'' (U '»»' : <0,„|'W|('l)|*»)+... '4.28! 



where 



S= lim <S(f,f ) (4.29) 



is the so called S-matrix. This is the basic result of time-dependent perturbation theory. 
For the case where 7i\ is time-independent we obtain for ((j) m j S |(j)„) 

<<|>m| S \tyn) 

= 5 m „ + <^) m |7Yi|^)y y df ie iaw ' 

df! / df 2 e''" (a, "- <0) e''' 2(a, - a, " ) +... 

= 8 mn -2ni8(E„-E m ) [<(j) m ] 7Yi |(j) n ) + <(j) m | 7iiGo(E n )7ii |<t>„> + ...] 
= S ffl „ - 27t/8(£'„ - £ m ) <(j) m I T + (£„) |0„) 



where we have used the relation 



The transition probability is then 



k4J<sk>l 2 = -Vi^m|wik>r / df 1 df 2e i<o ™"-'' <o ™' 2 + ... 



h 2 



\{§m\ 












|<<t>m 




|0n>| 2 



(4.30) 



y df exp ^ j = 2jt/?5(£) (4.31) 



f dfe lUhm ' J dt+... (4.32) 
df + ... 

so that the transition probability per unit time, W mn , is given by 

271 1 2 

W m „ = — <(j) m 7Yi (j)„) 8(E„-.E m ) + ... (4.33a) 
n 1 1 11 

= -r- (<|>m T + (E„) <|>„) 8(£„-£ m ) (4.33ft) 
n 

which is the well known Fermi's golden rule. Note that the higher order terms are included automatically by replacing 
Hiby r-(E n ). 
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Since |(j)„) is normalised, then |\(/ + (f)) is also normalised, hence 

1 = ( ¥ + (f)|y + (f)) = <4> B | A'(t, t a )A(t, t ) \ty n ) (4.34) 

and for the 5-matrix we have 

S f S = 55 + = 1 (4.35) 
which can be expressed in terms of the t-matrix using (4.30) 

<^|T + (£„)|^)-<^|T-(£„)|^> 

= -2tu^>«| |*m)<^| T~(E n ) \^ m )8(E m -E n ) (436) 

m 

This is equivalent to the optical theorem of scattering theory. 

It is convenient to rewrite the expression for S in the so-called interaction picture. 



S=1+ h 



-j j>-i-oo 

- J 6tin[(ti) 

(j) 2 J_ dhn[(h)J 1 dt 2 n[(t2) + ... (4.37a) 



%\(t) = e' H ° m ndt)e-' H ° m (437b) 
where we have used the relation 

gofa " f 2) = ®(f i " f2)g (f 1 " f2) 

= ®( h -t 2 )-e- m ^- h)lh (4 ' 38) 
n 

Note that the restrictions on the intermediate integrations in (4.37) can be dropped if the nth term is divided by n\. 
Also, in writing products of Wj's it is essential to preserve the original order of the operators. To make this ordering 
explicit we define the time-ordered product of operators 7i\ as 

T [H[(ti) ...n[(tj) ...] =H[(f l )H[(t 2 ) ...■H\(t n - l )7i\(t n ) (4.39) 

where t\ > t 2 > ... > f n -i > t n . Thus (4.37a) becomes 



S = J2 {jY^ I dfldf2 -dt n T[H\(.h)...H\{t n )] 



(4.40) 



which is just a compact way of writing the sum. 
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c) Application to Scattering Theory 

We consider Hq = p 2 /2m = -h 2 V 2 /2m, and Hi(r,r') = 8(r-r')V(r) is of finite extent (i.e. V(r) decays fast enough as 
r —> oo). Here Ti has a continuous spectrum from to +°°, but may develop discrete eigenvalues below zero. 

The problem of physical interest for E > is scattering. An incident article of energy E = h 2 k 2 /2m described by 
exp(/£ • r)/\^Q. is influenced by the perturbation V(r) and its wave function is modified. We must find the asymptotic 
behaviour as r — > °°. From (4.16) we have 

= -^=e ikr + -^= f dr! dr 2 G7 ) (r,r 1 )T ± (rur 2 )e ikr2 (4.41) 

Using the expression for the 3-D Gq obtained earlier we can write 

s/QyHr) = e ikr , / d 3 n d 3 r 2 7*(r 1 ,r 2 )e* r2 (4.42) 

4th 2 J \r~n\ 

where k = \f2mEI~h 2 . For large r we can write 
VQyHr) ^ e ik r 

tikr r 

— J d 3 n d 3 r 2 <r, 1 7\E) |r 2 > e* r > (4 43) 



m e ±ikr 
2th 2 ' 

iikr 



ikr m e 



= e - 



2th 2 r 



(±k f \ T'HE) \k) 



where kf is a vector of length k in the direction of r, and T = Q.T. 

In (4.43) \|/~ must be excluded as it describes a physically unacceptable incoming spherical wave. The quantity 
of importance is the scattering amplitude f(kf,k) which is defined by 



r— 



\|/(r) — > const. 



e* T +f(k f ,k) — 



(4.44a) 



ill 

f(k f ,k) = -—J (k f \ T + (E) \k) (AAAb) 
2th 

Thus the f-matrix is essentially the scattering amplitude, which is directly related to the scattering cross section 



do o 



m 



2 



47t 2 r 



(k f \T' + (E)\k) 



2 



(4.45) 



Substituting (4.44) into (4.1 1) we obtain 



/<*,*) = -^VOty-*) + / ^^^-f^k) (4.46) 

Thus to 1st order in V 

f(k f ,k)~-^V(k f -k) (4.47) 
2th 

where V(q) is the Fourier transform of V(r). This is the Born approximation for the scattering amplitude. 

Another derivation of (4.45): The probability per unit time for the transition & — ■> kf, W kf k, times the number of 
final states divided by the solid angle At and by the flux j = v/i2 of the incoming particle. 



da _ Q. 
di2 4tv 



J dE f N(E f )W kfk (4.48) 
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Using (4.33) we obtain the total cross section 



f da 



(4.49a) 



= -^£K*/l T W>l 25 ( £ /-^ 



v ft 

*/ 

= ^ E <*i T "(^ I*/) T+ (^ i*> § ( £ / - ^ 
= 1S [<*i t+ ( £ )i*>-<*i t "( £ ) i*>] 

= -^Im<*|T + (£)|&> 

= ^Im/(*,ifc) (4.4%) 

It 

which is the optical theorem which connects the total cross section with the forward scattering amplitude. 

As an example consider the case where V(r) is an attractive Coulomb potential V(r) = -e 2 /47t£or. Then / is given 

by 



/ = 



-tT(l - it) exp (lit In sin (±0)) 
1X1 + if) 2* sin 2 Q 9) 



k=,\^- lmfc>0 (4-50) 



t = 



me 2 



4n£(fi k 

where 8 is the angle between k and kf. The poles off occur when the argument of T(l - it) is a nonpositive integer, 
i.e. when \ — it = -p where p = 0, 1,2 or it= 1+p = n and n = 1,2, 3 Thus the discrete eigenenergies are 

which is the standard result (in SI units). If the potential is repulsive we must replace i by —i. in this case the argument 
in the F function in the numerator in (4.50) cannot become a non-positive integer; thus / has no poles, just as we 
would expect for a repulsive potential. 
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5. Tight-Binding Hamiltonians 

a) Introduction 

Consider the so-called tight-binding Hamiltonian of the form 

W = ^|/)e,</| + ^|Z>V, ffl <m| (5.1) 

/ lm 

where each state \l) is an atomic-like orbital centred at the site /. The sites {/} form a lattice. For the moment we shall 
consider the periodic case where the Hamiltonian is invariant under translation by any vector /. 
The eigenfunctions obey Block 's theorem 

l|/„*(r) = -j=e ikr u nk (r) (5.2a) 

Mf nk (r + l) = e ikl \f nk (r) (5.2b) 

where u„k(r) is a periodic function and k is restricted to a finite region of space, the Brillouin zone. 

It is useful to introduce Wannier functions associated with the band index n and the lattice site / defined as 

w n (r -1)=-L j2 e~ ik l Vnt(r) (5.3) 

vn k 

These functions form a complete set for all n and /. Any operator can be expressed in the Wannier representation. It 
may happen that the eigenenergies associated with a particular band index no are well separated from the rest. In this 
case the matrix elements of Ti between w„ a (r - 1) and w n (r - m), where n ± no, may be smaller than |e„ - E„ 1 . Then 
these small matrix elements may be omitted and the band no is decoupled from the rest. In what follows we work 
within this subspace and for simplicity omit the subscript no to obtain 

w(r-l) = {r\l) (5.4a) 

(l\n\m) = £,8 lm + V lm (5.4b) 

In the usual notation (with Vu = 0). Since Ti is periodic we can write 

Ei = £o (5.5a) 

V,„ = V/,- m (5.5b) 

We shall consider also the case of 2 interpenetrating sublattices where 

f £i if i e sublattice 1 . 

£i= „ ™u,-«:~o (5-6) 



£2 if i e sublattice 2. 

For simplicity we assume only nearest neighbour interactions 

y _ (V i j nearest neighbours ^ ^ 

1J 1 otherwise. 

Usually in the literature V is assumed positive. It is actually physically more realistic to consider V negative, since 
this guarantees that the energy k increases as E increases. The Green's functions for the 2 cases are related by 

G(l, m; E + is, {£,}, V) = -G(l,m; -E - is, {-£,}, -V) (5.8) 

From Bloch's Theorem we have 

By substituting (5.9) into Ti we can derive the eigenenergies 

^ XL e ' k l general case 

E(k) = £ + { 2Vcos ( fo ) 1_D (5 10) 

2V [cos(fcifl) + cos(fc2fl)] 2-D square lattice 

2 V J2d=i cos (^ fl ) 3-D simple cubic lattice 

In all cases -%/a < k < %/a which is the first Brillouin zone. Remember that the number of distinct ^-vectors in the 
Brillouin zone is equal to the number of unit cells in the crystal. Note that the bands go from -TV < E < TV, where Z 
is the number of nearest neighbours. Thus the band width is given by 2TV, but is often quoted as TV. 
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c) Green's functions 

The Green's function can be written 



Note in particular that 



, z-E(k) 

G(l,m;z) = ^-—j^- (5.11b) 

O f Jk(l-m) 

dk — (5.11c) 



N(2n) d J BZ z-E(k) 

Q. I dk 
AW* J BZ T^E{k) 



Gl/./:. > T77^ / —7^777 (5.\UI) 

JB2 



lim G(l,l;z) 



dEp(E) 



z->°° J z — E 

, Jp(E)dE (512) 



z 

1 



i) 1-D Case 



£ r+ic/a e ika(l-m) 

G(l, m; z) = tT ,„_^ / dk 
J-n 



N(2k) J_ K/a z-Eq- 2V cos(ka) 



— / diSf- 



2k J_ k z-eo-2Vcos(j) 

By symmetry the integral depends on the absolute value \l - m\. We transform it to an integral over the complex 
variable w along the unit circle 

-If w v ~ m] 
G(l,m;z)=—fdw (5.14) 
2mV J w l - 2xw + 1 

where 

and the 2 roots of w 2 - 2xw + 1 = are given by 



p± = x =F Vx 2 - 1 (5.16) 



where Vx 2 - 1 we denote the square root whose imaginary part has the same sign as Im(x). Note that p+P- = 1. Also 
|p + | < 1 and |p_j > 1 unless -1 < x < 1. Then both roots lie on the unit circle and the integral (5.14) is not well 
defined. This condition gives the continuous spectrum of Ti which lies on the real axis between £o - 2 V and Eo + 2V. 
By the method of residues we have 

-1 p + |; - m| 



yp+ - p r (5.i7) 

= p"- ml 

x /(z-£o) 2 -(ZV) 2 + 

For z coinciding with the spectrum we have 

+ Ti l i |/_m| 

G ± (lm;E)= - = (xT'Vl-x 2 (5.18) 

v/(Zy) 2 -(£-e ) 2 V J 
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The density of states per site is given by 



0(E) = T-lmGHUlE) 

0(ZV-|E-£o|) (5 ' 19) 



7tv/(ZV) 2 -(£-eo) 2 

Note the square root singularities at the band edges E = +ZV. Note also that G decays exponentially with \l - m\ when 
z does not coincide with the band, but does not decay when z is a band energy. 

ii) 2-D and 3-D cases 

In 2-D and 3-D there is no closed analytical solution. The integral corresponding to (5.13) above can be reduced 
to a form containing elliptic integrals. However there is an important feature which does not appear in 1-D, namely 
singularities within the band, so-called van Hove singularities. 

In 2-D the DOS has step-like singularities at the band edges and a logarithmic singularity in the middle of the 
band. 

In 3-D the DOS has square root singularities at the band edges and 2 step-like singularities within the band. 
A common approximate form for a 3-D Green's function is the Hubbard Green's function with the following 

form 

G(l, I; z) = , 2 = (5.20a) 

z-e + V(z-e ) 2 -(ZV) 2 

P(E) = ^^vW^M (5.20*) 
which has the correct band edge singularities but no internal features. 

iii) Bethe-Lattices and Cayley Trees 

A Bethe lattice or Cayley tree is a branching structure containing no rings. Strictly speaking the Cayley tree can have 
a variable coordination number, but we shall ignore this possibility. The method of solving for the Green's function 
on such a lattice is a useful introduction to a more general method. We characterise the tree by its branching index 
K = Z-\. K= \ corresponds to a 1-D lattice. 

Consider the tree to have been cut at some branch. Then the root node is connected to K branches and every 
other node is connected to K + 1. We denote by G'(z) the diagonal element of the Green's function on the root node. 
If we join K branches together with a new atom to form a new root node then the Green's function of the new node is 
related to those of the old nodes by Dyson's equation where Tio is the Hamiltonian for K trees plus a single atom and 
7i\ contains a matrix element V between the new atom and the old roots. Hence 

K 

G'(n, n; z) = G (n, n; z) + G (n, n; z) VG'(i, n; z) (5.21a) 

K 

= G (n, n; z) + G (n, n; z)V V G (i, i; z)VG'{n, n; z) 



!=1 



1-G («,«;z)V 2 J]Go(m;z) 



(5.21*) 

G Q (n,n;z) (5.21c) 



where i runs over the K root nodes of the old trees and Go refers to the Green's function before the perturbation V is 
switched on, and Go(«, «; z) refers to a single isolated atom with Go = [z - £o] -1 , where we can set £o = without loss 
of generality. When the trees are infinitely large we have G'(n, n; z) = Go(i, i; z), as the end of an infinite tree must be 
invariant. Substituting in (5.21c) gives 

G'(n, n;z)=[z- KV 2 G'(n, n; z)] ~ l (5.22a) 
KV 2 G'(n, n;zf- zG\n, n; z) + 1 = (5.22*) 
z±\fz 2 -AKV 2 

G (n, n; z) = (5.22c) 
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where the sign of the square root is chosen to make the imaginary part of G be the opposite of the imaginary part of 
z. In order to obtain the Green's function on a typical site of the lattice we must perform one more iteration which 
connects K + 1 rather than K trees together. Hence 

G(n,n;z)= [z- (K+ l)V 2 G'(n,n; z))~ l (5.23a) 



1 (*-l )z - T (* + 1)^4]^ {523b) 



4K (K+iyvt-z 2 



Note the peculiar property of the Bethe lattice: that the continuous part of the spectrum (where the square root is 
imaginary) is narrower than expected (-2^/KV < E < 2^/KV), but that there are discrete singularities at z = ±(K+ 1)V 
where the band edges would be expected. Note that the case K = 1 reduces to the 1-D case discussed above. 

The Bethe lattice is a useful model for amorphous semiconductors where it is possible to build in short range, 
chemical, order without an underlying crystal lattice. More recently it has been used as a boundary condition for cluster 
models of amorphous systems; the so-called cluster Bethe method. A more complex variation involves building a tree 
from rings of atoms such as rings of 6 4-fold coordinated atoms. This is termed a Husimi Cactus. 

c) The Renormalised Perturbation Expansion 

Consider again the tight-binding Hamiltonian Ti =Tio + Hi where 



n Q = J2\ l )Zi( l \ (5.24a) 
i 

ni = vY,V)(m\ (5.24Z?) 



lm 

where the sum in (5.24b) is over nearest-neighbour sites only. Then we can write 



G = Go + Gon + Q n iQ n iQo + ... (5.25a) 
G(l,m) = G (/,fw) + ^G (/,ni)<ni [Hi\ni)G (n 2 ,m) 

«1«2 

+ Go(/,fli)(Birfti|fl2>Go(fl2,fl3) 
n\ ...ri2 

X<«3 \Tii\n 4 )Go(n 4 ,m)+ ... 

(5.25b) 



From (5.24) it is clear that Go(«i,«2) = 8 ni „ 2 Go(ni) where 

G (n) = — !— (5.26) 
z-e« 

and (ti\ \H1\n2) is non-zero only when n\ and «2 are nearest neighbours. Hence (5.25) can be simplified to 

G(l,m) = 8 /m G (Z) + G (l)VG (m)8i, m±l + ^ G (/)VGo(«i)VGo(»i) + ... (5.27) 

We can keep track of all the terms in (5.27) by noting that it consists of a sum over all possible paths from / to m. Thus 
for a particular path we include 

i) a factor Go(/) for each site visited and 

ii) a factor V for each step from one nearest neighbour to the next. 

The most general path from / to m is found by decorating a skeleton path consisting of a self-avoiding path from 
/ to m. For each site ft i we can replace Go(«i) by the sum over all possible paths beginning and ending at iti, that is 
by G(n i,«i) Unfortunately, when we come to consider a second site «2 we must be careful to avoid double counting. 
Some of the paths in the decoration of have already been included in the decoration of ti\. Thus the decoration 
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of «2 consists of a sum over all paths from n 2 to n 2 which do not visit m. Let us write this sum as G(n 2 ,n 2 [iti]). 
The exclusion is included automatically if we set £„, = °° in the evaluation of G(n 2 ,n 2 [n 1 ]). Thus we can rewrite the 
summation in (5.27) as 

G(Z, m) = J2 G(l, l)VG(n u ni[l])VG(n 2 , n 2 [l, n^V ... 

(5.28) 

VG(m,m[l,ri\,n 2 , ...]) 

where the summation is over all self-avoiding paths starting from / and ending at m, Z — > n\ — > n 2 m. In 

particular for the diagonal matrix element we have 

G(l,l) = G (Z) + G(l,l)VG(n u ni[l])V ... G (l) (5.29) 

Note that the last factor is Gq(1) because all decorations of / have already been included as decorations of the initial 
site. 

It is worth noting at this point that in treating the Bethe Lattice we considered truncated trees. This is equivalent 
to the partial summations considered here. 
Equation (5.29) can be rewritten as 

G(Z, Z) = G (Z) + G(Z, Z)E(Z)G (Z) (5 .30a) 

where E(Z) is called the self-energy and is given by 

zZ(l) = J2VG(n u ni[l])V... V (530b) 

(5.30a) can be solved for G(Z, Z) to give 



1 



z-e,-Z(Z;z) 

It is sometimes useful to write down an expansion for E in a similar form to (5.28) or (5.29) 



(5.30*/) 



L (') =y E \, rm (5-31) 

^ z-e ni -Z(m[Z]) 

In the case of the Bethe Lattice E(«i[Z]) is simple the self-energy at the end of a truncated tree ie 

E(n 2 [Z > m]) = E(n 1 [Z]) (5.32) 

which can easily be solved to yield the same result as before. 

The self-energy has a number of useful properties which we may use later. Unlike the diagonal elements of 
the Green's function itself E is purely real when z is on the real axis close to a discrete eigenvalue of Ti. E has an 
imaginary part when z is real but in a continuous region of the spectrum. This also applies to localised and extended 
states respectively. Thus Im E provides a measure of whether the eigenstates can carry a current. 

Also, since everything except E in (5.31) is real (when Imz = 0) E can only acquire an imaginary part when the 
sum contains an infinite number of significant terms such that the sum does not converge in a conventional sense when 
all the terms are real. Thus in principle we can write down the complete series for (5.31) where all terms are real. 
However, the solution, as for example is implicitly implied by the self-consistent condition used to solve the Bethe 
Lattice, has a complex E. 

In a finite system the sum ends when all paths are exhausted and E remains real. In a localised system the terms 
involving very long and distant paths make negligibly small contributions. Again E remains real. 

A more mathematically precise way of formulating this condition involves using a complex z and studying the 
limit Imz — > 0. 

1' I E(Z { = ^ l° ca ^ sec * or discrete states 
inu->o± \ >0 for extended or continuous states 
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a) Single Impurity Scattering 

We now redivide our Hamiltonian into 2 parts where Tio describes an ordered system and Ti\ describes a single 
impurity at site /. Thus 

Tio = ^\m)Zo(m\ + V^\n){m\ (6.1a) 

m nm 

Wi = |/)e</| (6.1ft) 
We can study either the Greens 's function Q or the t-matrix T 

q = Qq+QqU x q q + g n + . . . (6.2a) 

T = Ti\+Ti\GoTix+TixGoTixGoTi\ + ... (62b) 

Substituting (6.1b) into (6.2b) gives 

T = \l)e </| + 1/> e </| Go \l) e </| + \l) e (l\ Go \l) e </| Go \l) e </| + . . . 
= |/)e (l + eG (/, /) + (eG (/, I)) 2 + . . .) (l\ 

^I^U)^ 1 {63a) 

G = Go + QqTQo 

Note that the poles of G (or T) coincide with those of Go with the addition of a pole at 

G (U;E P ) = - (6.4) 
e 

Note that this must be outside the band of Tio; otherwise Go has an imaginary part and (6.4) cannot be satisfied. 
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+ ... 



b) Two Impurities 

Now let us consider the case where 7i\ consists of the potential due to 2 impurities. That is 

7il = \l)z{l\ (6.5a) 

H m = \m)z'{m\ (6.5b) 

7iQi = 7iQ + 7ii (6.5c) 

T~iom = T~io + T~im (6.5d) 

with corresponding Green's functions Go, Goi, Go m , Q. Since we have already calculated Goi = Go + Go%Go we can 
consider 7Yo/ as the unperturbed part and 7i m as the perturbation. Thus 

G = Goi + GoiHmGoi + GoiTimGoiHmGoi + ■■■ 

e' (6.6) 

= Goi + Goi \m) — — (m\ Got 

1 - e Goi(m, m) 

By substituting mow for Goi we can derive the final result. However a more symmetric formulation can be more 
revealing. Let G' be the 2 x 2 submatrix of Go for the 2 sites / and m. Note that H \ is effectively a 2 x 2 diagonal matrix 
already. Then 

T = n 1 +n 1 G"H 1 +niG"H 1 G"H' 1 + ... 

e 0\ fe 0\ f G (l,l) G Q (l,m) \ fe 
e'J + \0 e'J \G Q (m,l) G Q (m,m) J \0 e' 

e OW G (U) G Q (l,m) We \ (6-7) 
e'J \G Q (m,l) G Q (m,m))\0 e' 

/ G (/,/) G Q (l,m) We 
\G Q (m,l) G Q (m,m) ) \0 e' 

Note that the terms of (6.7) consist of a series of steps between / and m interspersed with any number of steps which 
remain on the same site. This latter can be summarised by f; which is the t-matrix due to the single impurity at /. Thus 
each of the factors e or e' in (6.7) can be replaced by f; or t m and the G' lose their diagonal terms. Thus 

ti N ft, \ f G (l,m)\ ft! 
.0 t m ) \0 t m ) \G Q (m,l) J \0 t m 

't, 0\f G Q (l,m)\ft, 
t m J \G Q (m,l) J \0 t m 

G (l,m)\ ft, ' 
G (m,l) J \ t m 

We have thus separated the contributions to T into those due to scattering from a single site and those associated with 
multiple scattering between sites. 
The result can be written as 

T =fmi (% + T m + TiGoTm + T m G Q T,) (6.9a) 
where the common factor f m i is given by 

1 - t m tiG Q (m, l)G Q (l, m) 
where t m = z'l (1 - e'Go(m, m)) the value of the t-matrix for scattering from a single site. 

The various terms derived above can be given a physical interpretation. The Green's function G(l,m) represents 
propagation from site / to site m. Terms of the sort Gq(1, n)t„Go(n, m) represent propagation from I torn with scattering 
from site n, etc. Thus it is customary to call the Green's function G(l, m) a propagator. Sometimes the term locator is 
used for the diagonal elements, G(l, I), whereas others use the term propagator for all the matrix elements of G- 

To first order in t m we can write the t-matrix as 

T = T l + T m + 0(T t 2 , Tl 7{I m ) (6.10) 

and the multiple scattering terms can be ignored. In fact since they depend on Gq(1, m) which decays as \l-m\ — > °°, 
the single scattering approximation is good when the impurities are dilute (i.e. far apart). Note that the one exception 
is a 1-D system where multiple scattering cannot be ignored as Go(l,m) does not decay. 



T 



(6.8) 



+ ... 
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c) Infinite Number of Impurities 

In this section we consider the case where a finite percentage of the sites are occupied by impurities in a random 
way. An example would be a binary alloy A X B\- X , where a fraction x of the sites are occupied by atoms A with site 
energy £o + £a, and the rest are occupied by atoms B with site energy e<> + Eb- In general we don't know which sites 
are occupied by A or B, but only the probabilities. We are thus led to the concept of a random or disordered system, 
whose Hamiltonian is not known. If there are no correlations between the random variables {e'„}, their probability 
distribution is given by 

P ({£'„}) = ]Jp (e'„) (6.11a) 

n 

with 

p (e'„) = x8 (e'„ - e A ) + (1 - jc)8 (e'„ - e B ) (6. 1 lb) 

In a random system we are usually interested in the average (over all configurations) of the physical quantities. More 
specifically we should like to know the average of the Green's function {Q) defined as 



(G) = J d{e' n } P({e' n })G ({e'„}) (6.12) 
In general {§) cannot be calculated exactly. The Lloyd model constitutes one exception. When p(e'„) is a Lorentzian 

jc e„ 2 + r 2 

Then, because p(e'„) has 2 poles £'„ = ±iF, each integration over e'„ in (6.12) replaces e'„ by ±;T; thus finally 
{G ± (E)) = G Q (E±iT). 

Other commonly used distributions are the rectangular 



and the Gaussian 



P(e'„)<!-W for]e'„j<5W {6Ua) 
1 " otherwise 



! ( e? 



This Hamiltonian (often called the Anderson Hamiltonian) is the simplest one to model a random system. Even within 
the tight-binding approximation there are physically important aspects which it ignores. 

i) the off-diagonal matrix elements Vy may also be random variables (off-diagonal disorder). 

ii) the quantity e'„ may depend on other sites in its vicinity of n. 

iii) the random variables {£'„} may be statistically correlated (short-range order). 

iv) in amorphous systems such as glasses the disorder is a associated with the absence of a lattice rather than the 
values of the matrix elements (topological disorder). Such disorder is much more difficult to treat. N.B. Some 
authors (e.g. N.Rivier) reserve the term topological disorder for systems containing odd numbered rings of atoms. 
This is often, but not always, the case in real glasses. 

In spite of the approximations built into the models it is still necessary to make further approximations to solve 
these models for {§). 
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d) Virtual Crystal Approximation (VCA) 

We use again Ti = Tio + Ti i where Tio is periodic and Ti i contains the random part. Then 

(G) = Go + Go(H l G) (6.15) 
since Go does not depend on the random variables. We assume 

(niG)~(ni)(G) (6.16) 

Taking into account that 

(n l ) = Y,\™)(e'm)(m\ = (e' m ) = e (6.17) 

m 

since (e'„) is independent of m. 
Substituting in (6.15) we have 

<Q(E)) = Go + GoZ(Q(E)) 
1 

" E-n -E (6-18) 
1 



E-(H) 

This approximation for (Q) consists in calculating (G(Ti)) as G{{Ti)) and simply shifts the energy levels by e = (e'„). 
This shift can always be made equal to zero by Tio = (Ti) and Ti\=Ti- (Ti). From now on we adopt this convention. 

The VCA fails completely for large values of (e„). However it seems to work surprisingly well for some alloys 
(e.g. Al x Gai- x As). 

e) Average t-Matrix Approximation (ATA) 

If T is the t-matrix associated with Tio and Ti \ we have 

{G) = Go + G Q (T)Go (6.19) 
where T is a complicated function of the individual f m 's 

T=/({f m }) (6.20a) 

t m = t4?7 T ( 6 - 20 ^) 

1 - e m Go(m,m) 

The basic approximation of the ATA is to put 

CO~/({<U}) (6.21) 

Because (6.21) has the same form as (6.20a) with an jw-independent (t m ), it follows that (T) is the same as the t-matrix 
associated with Tio and a periodic Ti \ = E, where £ is chosen to satisfy 

{tm) = : V „ E , r (6.22) 

1 -L.Go{m,m) 

Thus the energy is shifted by the complex self-energy E to give 

(G(E)) = Go(E-Z) (6.23a) 

where 

1 +{t m )Go(m,m) 

In the limit of weak scattering € m — > 0, we obtain, in lowest order that 

Y. ^ Ul) Go{m,m) (6.24) 



which is the correct limit. The ATA is also successful in the very dilute limit, where the concentration of impurities is 
< 1. 
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f) Coherent Potential Approximation (CPA) 

As in the ATA case, the CPA calculates (G) through an effective Hamiltonian Ti e , which in the simplest case is 
characterised by a single, complex, energy-dependent self-energy E. The CPA is based upon expanding Q in terms of 
Q e = (E- n e T l and U\ = U - 7i e . Thus 

G = Ge + GeTQ e (6.25fl) 



where the t-matrix T is 
with 



T'=f({f m }) (6.25b) 
(e' m -E) 



1 - (€ m - E) G e (m, m) 

As for the ATA the central approximation is 



(6.25c) 



(T')~f({(t' m )}) (6.26) 

Taking the average of (6.25a) and using (G) = Q e we obtain {T) = 0, which in view of (6.26) implies (t' m ) = 0. From 
(6.25c) we can write an equation for the self energy in the form 



' i u",„ ek;,.,,,,,,,)/ (6 - 27 " ) 

1 = ( 7 ) (6.27b) 



1 -(e m -E) G e (m,m) 

In (6.27) the quantity G e (m, m) depends on E through 

G e (m,m;E) = G (m,m;E-T.) (6.28) 

The CPA has been impressively successful. In the weak scattering limit it reduces to (6.24). It behaves correctly in 
the strong scattering or dilute limit, and interpolates properly between them. 

It does however have failures. The only approximation used is that (t' m ) = implies (T') = 0. The physical 
meaning of this can be found by expressing T in terms of the T' m 's as 

f = T 'm + Y, T ^ T » + £ ^ ' mG' ()T ' mG' ()T ' r +... (6.29) 

When the {e' m } are independent random variables the quantities {t' m } are also independent random variables. Thus 
the 1st 3 terms in (6.29) depend on (t' m ) alone, the 4th term contains contributions like 

(6.30) 



Yl 1"^ ( f ») Go(n,m)Go(m,n)Go(n,m)(m\ 



Similarly higher-order terms will involve terms in (tj£) for any integer n, which correspond to multiple scattering 
from a fixed number of sites. We conclude then that CPA incorrectly treats such multiple scattering. Such scattering 
is of great importance for eigenstates which significantly enhanced in the vicinity of a particular cluster. The multiple 
scattering is what traps the electron in the region for a long time. Hence at energies where the DOS is dominated by 
resonance or localised states, the CPA is expected to fail. Since the tails of the band are composed of such states the 
CPA is a very bad approximation there, in that it predicts no tail at all. The CPA tends to eliminate the structure in the 
DOS associated with such states. If the probability of occurrence of a special cluster is small the CPA is good. Such 
special clusters will tend to be less likely in 3-d than in 1-d as the number of atoms involved will be larger. Thus 
the CPA is better the higher the dimensionality. In this it has much in common with a mean field theory of a phase 
transition. 
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g) Extensions of the CPA 

In our derivation of the CPA we assumed that the random part Ti\ was a sum of local terms. This was a necessary 
assumption. It can however be easily generalised to the case where Ti i involves a finite number of local orbitals. In 
this case we can write 

n x = £^|m,v>e' JB ,w'<v',m| (6.31) 

m vv' 

and the effective Hamiltonian takes the form 

m vv' 

where the matrix E V v' is determined by a matrix equation of the form (6.23b) for the ATA and (t' m ) = for the CPA. 
Care must be taken with the order of the matrices. 

Another example is a special case of off-diagonal disorder where 

V mn = V m + V n (6.33) 

then Ti i is the sum of local terms of the form 

|m>^(e' ffl 8 ffl „ + V ffl )<«| (6.34) 

n 

where the summation is over the nearest neighbours of m. 

This is the case for the special binary alloy system where 

V AB =\(Vaa + V bb ) (6.35) 
and lattice vibrational systems with random spring constants where 

t m =-^V mn (6.36) 

n 

When (6.36) is satisfied, Ti i can be decomposed into bond contributions of the form 

\n) V nm (m\ + \m) V mn (n\ + \m) V mn (m\ + \n) V nm (n\ (6.37) 

The CPA can also be generalised to other representations than tight-binding. For example, the so-called muffin-tin 
CPA avoids the problems associated with off-diagonal randomness in the tight-binding case. 

These generalisations of the simple CPA share with it its advantages and difficulties. The limiting cases are 
correctly reproduced and the resulting DOS has proper analytical behaviour, i.e. no negative DOS. However this 
was only achieved by assuming special forms for the Hamiltonian. There have been numerous attempts to develop 
cluster-CPA methods. These involve for example: 

a) Forcing the system to have the periodicity of a larger supercell rather than that of a single atom per unit cell. 

b) Embedding the disordered system in an ordered effective medium. This category includes certain versions of the 
Recursion Method and the Method of Moments as well as the Cluster Bethe method. 

7. Electrical Conductivity 

Disorder has a much more profound effect on transport properties than on the DOS. The DC electrical conductivity is 
no longer infinite (as T —> °°), and we can even find a metal-insulator transition. 
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a) Definition and Basic Results 

An electric field E induces a current density j. The conductivity is defined as the coefficient of the linear (in E) part of 
j- 

j a (r,t) = j dx j 6r'a a? (.r,r';x)E^r',t-x) (7.1) 

where the subscripts a, p denote Cartesian coordinates and a summation is implied over the repeated index p. We shall 
assume that E and j are both along the x-axis, so that we only need to consider a**; for simplicity then we drop the 
subscripts. 

Usually E and j vary slowly over distances of the order of Iq, where Iq is determined by the condition a ~ for 
jr - r'j >> /(). In this case we can perform the integration over r and the average over r to obtain 

j(t)= / dxa(T)Ea-T) (7.2) 
Jo 

where 

a(x) = ^ J dr dr' a(r, r'; x) (7.3) 
We now consider o(x) or its Fourier transform 

o(co)= / dxo(x)e iayc (7.4) 
Jo 

If £(f) is given by 

E(t)=Fe- i(at + F*e i(at (7.5) 

we have 

j(t) = c((a)Fe- i(0t + o(-(j))F*e i(0t (7.6) 

the reality of j requires that o(-co) = 0*(co) from which it follows that the real (imaginary) part Oi (02) of o is an even 
(odd) function of co. because of causality, o(x) is nonzero only for x > 0. Consequently, o(co) is analytic for Im co > 
with the possible exception of co = 0. Taking this into account as well as the fact that o(co) — > as co — > 00 (see later), 
we can see that 

^'^=0 (7.7) 



CO +W-C0 



To prove this consider a contour around the positive half plain. There are no poles. By taking the limit s — > + and 
by considering the imaginary (real) part of (7.7) we obtain 

1 f dco'cWco') 

Oi(co)=-P/ - (7.8a) 

n J co - co 

1 / dco'oi(co') A 

o 2 co = — P / + - (7M) 

% J CO -co CO 

where iA is the residue (if any) of o(co) at co = 0. Thus, knowing Oi(co) for co + one can calculate 02(0) and 
consequently o(co). The constant A can be obtained from the behaviour of o(co) at infinity. 
(7.8) are usually referred to as Kramers-Kronig relationships. 

The simplest and crudest way to obtain o(co) is by using Newton's equations for the electronic drift velocity v in 
the presence of the field Fexp(-/cof) and a friction force -mvli u (x tt is the transport relaxation time): 

mv 

-mmv = + eF (7.9) 

X tt 

The current is given by j = nev where n is the electronic density, we obtain for o(co) 

a(«»~ , (7.10) 

m(l - (COXtj-) 
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For co — > °° the electronic motion is classical, all scattering is negligible, and only the electronic inertia ( as measured 
by its mass m) matters. Hence 

2 

oo->~ ne 

a(co) — > i (7.11) 

mCO 

A slight generalisation of this approach, capable of treating non-spherical Fermi surfaces, consists in considering the 
whole equilibrium distribution transposed rigidly in &-space by an amount 8k = mv/Ti, where v is again given by (7.9). 
A more sophisticated approach is one given by Boltzmann's equation, which determines the electron distribution in 
£-space in the presence of an applied electric field and scattering mechanism, which tends to restore equilibrium. The 
advantage of this approach is that the phenomenological parameter x tt is related to the scattering potential V(r) by 

— = ^pF»im P ^ / Sdq V\q)(\ - cos 9) (7.12) 
x tt h y 4% J F 

where Pf is the DOS per unit volume per spin direction, n lmv is the number of scatterers per unit volume, the integration 
is over all directions of q (E(q) = Ef) and V(q) is the Fourier transform of V(r). (7.12) follows from Fermi's golden 
rule for the transition rate from a given state on the Fermi-surface to all others. The extra factor (1 - cos 8) accounts 
for the fact that what matters for transport is the momentum change along the direction of initial propagation. If V 2 (q) 
is isotropic ( or more generally , if it does not contain a /^-spherical harmonic) then the cos term averages to zero, 
and x tt = x, where x is given by (7.12) without the cos 9 term. 

b) General Formula for the Conductivity 
i) Linear Response Theory 

Consider the density matrix p(r,r') which is defined as 

p(r,r') = Y,Ur)V a (r') (7.13a) 

Otocc 

= ]T<r|a><a|r'> (7.13ft) 

Otocc 

where the sum is over all occupied states. The density matrix is closely related to the Green's functions and shares many 
of their properties. A more general definition replaces J2a xc by J2af<*' where f a represents the value of Fermi-Dirac 
distribution function for the state |a). 

We shall study the evolution of the density matrix using the equation 

-ihp=[H,p] (7.14) 

We write the unperturbed equilibrium Ti and p as Ho and po respectively, and the corresponding deviations as 8Ti(t) 
and 8p. We assume further that Jio and po are time-independent, whereas bTi(t) is explicitly time-dependent. Then 
to 1st order we can write 

-ifi (p + Sp) = [H Q , po] + Po] + [H Q , 8p] (7.15a) 

-ihtp = [8H(f), po] + [Ho, 8p] (7.15ft) 

By taking the matrix element of both sides with eigenstates |a) and | p) of Tio we obtain 

-iti (a |5p| p) = (a \ [&H(t), po] | P) + (a | [Ho. §p] | P) (7.16a) 
= -(/a-/p) (a|5W«|P) + (e a -ep) (a|5p|p) 

(7.16ft) 

which can be rewritten in the form 

| (a |Sp| p) - l - (e a -e p ) (a |8p| p) = -- (f a -fy) (a \dn(t)\ p) (7.17) 
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Using the integrating factor method the solution of (7.17) can be written down 

where 8co a p = (e a - e p )/7j and 8/„p =f a -/p. 

We now specialise to the case where STi(t) = &H exp(r|f + mt). The rationale behind this is that the perturbation 
should be switched on infinitesimally slowly at t = -°° and we are interested only in long term rather than transient 
behaviour. We shall later take the limit r| — ■> + . Substituting this form into (7.18) and integrating gives 

a Sp p) = --8/op X '. J 1 (7.19) 
h p co-;r|-8co a p 

This is the basic result of linear response theory. 

In order to calculate the conductivity we require 

e (r\+i<a)t 

m(t) = -E p (7.20a) 

= ex ■ Ee^ +ia)t (7.20b) 

= L(-P8p) (7.20c) 



(7.20J) 



= "7^7 Y |<« m P>| 2 7^ 5 T Ee^' 

Llm z h 1 1 (r| + /co) (co - ir\ - 8co a pJ 



The real part of the conductivity o is found by removing the time-dependent electric field factor from (7.20d), taking 
the limit r| — > + , and expressing the result in terms of delta-functions to obtain 

^(«) = -^EK^lP)| 2 |^S(co-5co a p) (7.21) 

ap 1 

Note that this function is always positive: the ratio 8//8co is always negative. 
An alternative derivation of (7.21) is given by Economou. 
This formula is usually referred to as the Kubo-Greenwood formula. 

Using the Kramers-Kronig relations (7.8) we can now derive the expression for 02(a)). Combining these 
expressions we obtain finally 

I e 2 n g2 1 / 1 - 1 d\ I 2 ^faB I ' 



m Hm 22 -— ' lv 1 x| 71 h&(£> a R j (£> + is 

/ (7.22) 



ap 

e 



2 - 5/ a p 



' Qm 1 \ p x\P/\ ^8co aB co + is + 8co aB 



ap 

The term ie 2 nlm(0 = a</(co) is the so-called diamagnetic contribution to o(co), and the rest, O p (co), is the paramagnetic 
one, which can be recast as 

a„«D) = +i 7 r^- T I (« 1^1 P) I 2 (7-23) 

7 Q.m 2 h(s) ^ 1 x CO + w + 8co aB 

ap K 

The local conductivity a(r,r'; co) can be obtained by replacing the velocity operator matrix elements, p x /m, by the 
current operator matrix elements: 
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| (« \Px\ P) 1 2 l&^K$(r)K$( r ')> where 



-ih 



3\|/p(r) ,dy„(r) 



3x 



8x 



(7.24) 



where \|/ a (r) = (r|a). To obtain we must replace | (a \p x \ P) | 2 by 



a 



P) («WP). 

A useful identity relates the position operator x and the momentum operator p x : 



{a\p x \^) = l ^{a\YH,xW 



= /m8co a|3 <a|f| p) 



(7.25) 



It should be noted that this relationship generally fails in the presence of periodic boundary conditions where the 
position operator is not well defined. Nevertheless when all states are localised in a region much smaller than the 
period (7.25) is often used. 

Mott and Davis give an instructive way of deriving (7. 10) ( for co = 0) starting from (7.21). Their main assumption 
is that for weakly disordered systems the eigenfunctions have essentially constant amplitude (as in perfectly ordered 
systems) but phase coherence is maintained only over a distance of the order of the transport mean free path l u = v F i u . 
If the total volume D. is divided into cells, each of volume u, = lf r , then 



11. JVi 



(7.26) 



Multiplying (7.26) by its complex conjugate and cancelling the cross terms due to their random phases, one obtains 

2 



K«|A|P)| 2 = Z] / ^VaPxV 
i \Jy>i 

= -\J df-VoPxVp 



(7.27) 



Within the volume V, \|/ a and V|/p are assumed to be plane waves so that the integral can be performed explicitly. 
Substituting the result in (7.21) one obtains 

e 2 nl u 



2mvp 



(7.28) 



which coincides with (7.10) except for a factor of 2 (for (0 = 0). 

Thouless has also obtained (7.10) starting from (7.21) by assuming that \|/ a and are linear combinations of 
plane waves with coefficients which are uncorrected Gaussian random variables whose variance is 



l ti k 2 Q. 



(k-k a ) 2 + l 2 ti 



(7.29) 
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c) The Conductivity in terms of Green's Functions 

There are several different formulations of the conductivity in terms of Green's functions. Firstly we express o(co) 
in terms of the causal Green's function, g(E), which is important in many-body theory. For the present purpose, 
non-interacting electrons, g(E) can be defined as G~(E) for E <E F and G + (E) for E > E F , i.e. 

g{E) = [E+ is€(E - E F ) - %Y l (7.30) 

At T = we have 



§/a|3 



— = -1/ dE(a\g(E)\a)^\g(E + h( S )W (7.31) 

+ IS ZKl 



e a p + ho) 
Substituting in (7.23) leads to 

ifp'Yl f (\E 

0(w) = ^ + ^y ^{pjm^m} (7.32) 

where the Tr operation includes a factor of 2 for the spin. By using standard many-body techniques (7.32) can be 
generalised to finite temperatures. 

It is probably simpler to express Oi(co) in terms of Q = Q + - Q~ . for that purpose one notices first that 



-Sco, 

Substituting (7.33) in (7.21) we have 



2-8 (^to + 8e„R) = ( &E 8(£ - E a )8(E - E* + ha /— f{E + m (7.33) 
taB J to 



Dm 1 nco 

x Tr {p x 8(E + Tm- %)p x h{E - %) } 
e 2 K f + °° f(E)-f(E + h<£>) 



r dE f(E)-f(E+i 

J-00 Ti& 



(7.34) 



Tldm 2 

x Tr { p x \mQ + (E + Tm)p x Im Q + (E) } 

where again the Tr includes a factor of 2 for spins. Oi(co) can be expressed in terms of Q ± by noticing that 
\mQ + = i\{Q--G + ). 

Note that in the limit co — > these equations can be simplified using 

,. /(£)-/(£ + to) df 

hm h-> - (7.35) 

(0^0 h(£> 6E 

and that in the limit T — > this further simplifies to 

-^8(E-E F ) (7.36) 

Another useful form is found by transforming the p x operators to x operators using the relationship ihp x = m [H, x] . 
Substituting this in (7.34) and rearranging gives 

2e 2 

01 (0) = — Urn Tr {r| 2 (Q(E + iy\)xQ{E - ir\)x) 
{Q{E+iy\)-Q{E-iy\))x 2 } 

This latter form has proven particularly useful as the starting point for the development of a numerical method for the 
calculation of the conductivity. 
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In all this it is necessary to average over the ensemble of disordered systems, usually employing the arithmetic 
mean of o. There is much evidence that the distribution of In o is better behaved (approximately log-normal) and it is 
often argued that exp (In o) is more representative of the whole distribution of o's than (a). 

In order to calculate the arithmetic mean it is necessary to calculate the average of a product of 2 Q's, which is not 
generally equal to the average of the product. However we shall consider the results of such a calculation. Working in 
the ^-representation where both p and (G) are diagonal and considering the d.c. conductivity at T = we have 





Px 






m 





\\mQ + (E,k)\ 



(7.38) 



Taking into account that (k ](/? x /m)| = v x (k), and that 

G + (E,k) = [E - E(£) - E(k)]~ , where E = Ei + zE 2 is the self-energy we obtain 



2e 2 Ti v-- 



E 2 
^2 



~k ([Ef- Ei -E(k)] 2 + E 2 ) 
The summation over k is easier if we introduce an integration over 8(£" - E(k)) dE', so that 



(7.39) 



o(0) = 



2e 2 h 



dE' 



Y,v 2 x (k)W-E(k)) 



^2 



-£(*)] 2 + 



(7.40) 



The quantity in brackets, [], only depends on the form of E(£). For cubic symmetry lattices v 2 . can be replaced by v 2 ld. 
Then the k summation can be replaced by an integration over the Fermi-surface, so that 



\v(k)\ 2 



1 



(7.41) 



(2n) d dh 



5(£')v(£') 



where S is the are of the surface of constant energy E" and v(£") is the average of the magnitude of the velocity over 
this surface. 

In the weak scattering limit , where E2 is small , and Sv can be taken as a constant, the integration over E' can be 
performed explicitly to give 

e 2 VF Sf 



0(0): 



Pal (2%) d d 



For a spherical Fermi-surface this becomes 



1 Vf 1 
c(0) = <? 2 n- 



fe F 2 Pal 



(7.42) 



(7.42) 



where n is the density of electrons at the Fermi-level. Comparing this result with the classical one (7.10) we can 
identify x with h/2 |E 2 |- 

Improved results can be obtained by including so-called vertex corrections for (GG) - (G) 2 - However we have 
already seen that the CPA omits much of the important physics of disordered systems. 
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d) CPA and Vertex Corrections 

Consider a quantity of the form (QBQ). In the framework of the CPA this can be written as 

(G&G) = GeBGe +{(G-Ge)B(g- Qe)) 

= (GeTVeBGeTVe) (7-43) 
= (GeTGe) 

where the residual t-matrix is non-diagonal in the site indices. We now use the property that 

T'(m, n) = \m) f m 8 mn (n\ + ^ \m) t' m {m\ G e (m, m')T(m, n) (7.44) 

m'*m 

and substituting this in (7.43) to obtain 



T m „=(f m (l+J2 G e (m,m')T'(m',m") J G e BG e 

\ \ m'4m / 

x ( 1 + r in",n')G e (n',h) j t' n \ 



(7.45) 



The CPA for the vertex part T mn is to ignore the correlations between the f' m 's and the T'(m,m'ys in (7.45). Thus we 
can write 



( t'm ( 1 + X] G e (m,m')T(m',m") GeBGe 

\ \ \ in ' -f-m / 

x ^l + ^ r(n",n')G e (n',«)j^f'«^ 



(7.46) 



Since the t' m are independent random variables with zero mean (due to the CPA condition) the RHS of (7.46) is 
non-zero only when m = n. Hence 

r m „ = 8 m „r„ = 8 mn |n) y„ (n\ (7.47) 

Of the 4 terms inside the average in (7.46) 2 are proportional to (T) and hence are zero within the CPA; the 3rd is 
simply GeBGe, and the 4th is 

J2 Ge (T\l l')GeBG e T(ri, »)) Q e = G e(™< n)T n G e {n, m) (7.48) 

l,nlm ntm 

Usin S J2„^ m r « = J2„ r n~ T ffl we can rewrite (7.46) as 

J„ [l+{f„f n )Ge(n,n)Ge(n,n)} 

= Y, < f JJ Ge{nJ) i B "» + ^ 5 «J Ge{m ' n) (7 ' 49) 

lm 

In order to simplify this expression we define the quantity u(z, z) as 

«(z. z) = (t' n (z)f n (z)) [1 + (t'„(z)f„(z)) Ge(n,n;z)Ge(n,n;z')] (7.50) 
so that (7.49) now becomes 

Jn(z,z) = u(z,z)Ge(n,l;z) [B, m +y,(z,z')Si m ] G e (m,n;z') (7.51) 

lm 
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The solution of (7.51) can be found by making the substitution 

In = JnmG e (ml)Bi k G e (km) (7.52) 



mlk 

so that (7.51) can be rewritten in the form 



- uh nm - ^2 uG e (n,l)G e (l, «)Y/„, 



(m\g e Bg e \m)=0 (7.53) 



The only B-independent solution is when the quantity in [] brackets is itself zero; that is 

J2"G e (n,l)G e (l,n)ji m (7.54) 

The quantity y nm corresponds to the vertex-correction in many-body theory. (7.54) can be solved in k space to yield 

y(q;z,z) = u(z,z') [l - u(z, z)A(q; z, z)]~* (7.55a) 
Mq; z, z) = e iqm G e {0,m; z)G e {m, 0; z) 



k 



J nm (z, Z) = ^Y, e iq< "~ m) y(<i; z, z) (7.55c) 



N 

9 

These results can also be written in a diagrammatic representation. Note that for weak scattering mi— > (e n 2 ) and 

Ymn 1 > M&mn- 

If we are interested in the conductivity tensor we first note that B = p and that p in the £-rpresentation is an odd 
function of k due to time-reversal symmetry. Hence from (7.52) y„ = and the CPA vertex corrections vanish. Note 
however that for quantities like 

(G(m,n;z)G(n,m;z')) = -J- V e~ ik ^ - (7.56) 

N ^ l-u(z,z)A(k;z,z) 

k 

which is related to the transmission coefficient, they constitute a major contribution, 
e) Vertex Corrections beyond the CPA 

Clearly we have omitted some important terms from the CPA; (7.46) is an approximation. The most important of 
these form a group very similar to the CPA ones, and may, therefore, be equally important. They are best expressed 
diagramatically using the so-called totally crossed diagrams. In conventional notation we have 

Tr {A (g(z)BG(z)) } - Tr {A (g(z)) B (g(z'))} 

(7.57) 

= 2^ A [s G e (s,m; z)i mn G e (n, i; z)B ij G e {j,m; z)G e (n, I; z) 

and the equation corresponding to (7.54) is expressed in expanded form as 

Ynm = G e {n,m)uh mm G e {n ,m')uh nn ' 

+ G e {n, l)u8 nn G e (l, m)ub lv G e {n, l')uh mm G e {l ', m') + ... (7 58) 

IV 

This is really the same as before except that half the Green's functions have been turned round. That is G e (m, n)G e (n,m)^G e (m, n)G e (, 
We note that in the presence of time-reversal symmetry G e (m, n) = G e (n,m) but that a magnetic field will destroy that 
symmetry. Hence Ymn = Jmn- The difference between the two cases lies in the way the vertex part y mn is connected to 
the rest of the system. In the CPA case the 2 G e 's which connect y mn to A and B start and end on the same site, whereas 
in the case of the totally crossed diagrams the sites are different. The symmetry which made the CPA correction to 
the conductivity zero is no longer effective and a finite contribution results. This contribution plays a crucial role in 
localisation (see below). 
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8. Disorder and Localisation 

Earlier we tried to calculate the conductivity of a disordered system by assuming that the amplitude of the wave 
function is constant and the phase is random. This is quite a common approximation in the metallic regime where 
the phase may be considered random on length scales larger than the mean free path. However there is no reason in 
principle why the amplitude should not be affected. 

Consider the simple case of 2 coupled pendula. The transfer of motion from one to the other is facilitated by 
strong coupling and opposed by frequency mismatch. If we now consider an array of such pendula, then it is not 
difficult to imagine a region of strong frequency mismatch which a wave propagating in such a medium will find 
almost inaccessible. In fact one could imagine a wave surrounded by such regions which is confined to its own region, 
unable to escape. 

Previously we considered a single or double impurity. These gave rise to discrete states outside the continuous 
spectrum of the rest of the system. What happens then when we go over to a finite density of impurities? We should 
expect the discrete spectrum to be smeared out to form a continuum which forms a tail of the main band. The question 
arises whether such states remain localised or become extended. 

The answer to this can be found in the framework of the random tight-binding model often termed the Anderson 
Model. It contains in its simplest form the essential competition between the coupling V and the energy (or frequency) 
mismatch, which is characterised by a width 8e (8e = \£a -Eb\,T,W for binary, Lorentzian and rectangular distributions 
respectively) of the probability distribution of the random variable e„ . Thus the important parameter for localisation is 

Q=l (8.1) 

As we shall see an equally important parameter is the dimensionality d. Other aspects such as the shape of the 
probability distribution and the type of lattice are of secondary importance and are believed not to influence the alleged 
universal features of the problem. 

We now summarise the main results in the field. These results have usually not been rigorously proven but are 
generally accepted. 

i) There is a critical dimensionality d = 2. For d < 2 all eigenstates are localised, no matter how weak the disorder 
is. For d > 2 and for weak disorder the tails of the band consist of localised states whereas the interior corresponds to 
extended states. These regions are separated by critical energies, called mobility edges. As the disorder is increased the 
mobility edge tends towards the centre of the band until eventually, at the Anderson Transition, all states are localised. 
This can be represented on the E - W plain, by a mobility edge trajectory which separates the regions of localised and 
extended states. 

ii) For d < 2 the vertex correction to the conductivity, for weak disorder, is given by 

JULm 1 

where L M is dominated by the shortest of several upper cutoff lengths which may be present in the system, and L m , the 
lower cutoff length, is believed to be very close to the mean free path l u \ 

L m ~y/D^ = yVd (8.3) 

where d is the dimensionality, D is the diffusion coefficient which is related to the conductivity by the Einstein relation 

a = 2e 2 p F D (8.4) 

where p F is the DOS at L F per spin per unit volume (area or length). 

Examples of upper cutoff lengths L M , besides the geometrical dimension L, are; the diffusion length during the 
inelastic relaxation time v _ 

U = v^; (8-5) 

the diffusion length during the time CO -1 , where co is the frequency of an external AC field 

La, = VDh£>. (8.6) 
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The presence of an external magnetic field B introduces 2 characteristic lengths: 

L' B = s/DI^ (8-7) 

where ©o = eB/m, the cyclotron frequency, and the cyclotron radius 

L H = (8.8) 

For weak disorder, Z/ B L H , and hence L H is the relevant length. If we assume that the various rates (i.e. 1/x) are 
independent of one another, the effective upper cutoff length due to the combination of several is 

L-^ = L- 2 + U T 2 + L-2 +... (8.9) 



Ford = 2 (8.2) yields 



while for d = 1 (where L m can be taken as zero) 



So, =-4 m ^ (8.10) 



8o!~-L M (8.11) 



Equations (8. 10,8. 1 1) show that the corrections do not disappear as the upper cutoff length L M becomes very large. In 
fact they increase with L M until O — > 0. Thus truly metallic behaviour is not possible for d < 2. This is, of course, 
consistent with the statement that all eigenstates are localised. This fact is often masked by the fact that at high enough 
temperatures L T is quite short. Additionally the localisation length can be macroscopic for quite significant disorder 
and astronomical (literally) for weak disorder. 10 10 has been published. 

iii) The vertex correction (8.2) comes from the post-CPA terms discussed earlier. The CPA terms although quite 
similar produce no correction to the conductivity. The presence of electron-electron interactions eliminates this 
cancellation. In fact they produce a very similar result to (8.2). The cutoff L M now depends on the frequency co and 
the temperature T as 

L^~k B T/hD + (a/D (8.12) 

However in this case Lm does not depend significantly on the magnetic field. This is because the CPA terms involve 
Gnn£*mn which is unaffected by the field, while the post-CPA terms contain G* m G~ m which acquires an extra phase 
factor in the presence of B. 

iv) The theoretical results presented in (ii) and (iii) are in impressive agreement with experimental data, especially 
for2-d. 

v) It has been suggested that (8.2) can be used even for 3-d systems as we approach the mobility edge from the 
metallic side. For d = 3 (8.2) gives 

So, ~ 1/L M (8.13) 
When L M is a magnetic length this predicts a negative magnetoresistance proportional to B m . 
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a) One-Dimensional-Case 

Mott & Twose were the first to propose that all eigenstates of a 1-d system are localised. If we assume that each 
back-scattered wave is lost due to destructive interference, then this picture suggests that the localisation length X, 
defined by the relation 

A," 1 = - lim ({in |\|f m A|/ |) lm) (8.14) 

is about the same as the mean free path I; \\i m is the amplitude of the eigenfunction at m. 

Borland showed that for a random 8-function array the solution of the differential equation with fixed boundary 
conditions at one end grows exponentially on average with the distance. This exponential growth was seen to be the 
result of phase incoherence. The picture that emerges on the basis of this proof is that at every energy there are 2 
independent solutions growing exponentially in opposite directions. At particular energies (the eigenvalues) the left 
growing solution matches up with the right growing solution to form a localised eigenstate. According to this picture 
the localisation length is the same as the rate of growth of the solution of the differential equation. 

The exponential growth can be shown rigorously using the transfer matrix technique. The transfer matrix is a 
2x2 matrix which connects e.g. \\t, d\\tldx at x with d\\ildx at x . The concept of the transfer matrix transforms the 
propagation of the wave to a product of random matrices. There are exact theorem's which show under quite general 
conditions that the product of random matrices grows exponentially. (See e.g. various papers by J.B. Pendry and 
P.D.Kirkman). 

The transfer matrix technique can be used for the calculation of the transmission (|f| 2 ) or reflection (|r| 2 ) coefficient 
associated with a disordered segment of length L. Using Furstenberg's theorem it is easy to prove that 

1 <ln|f| 2 > =-| (8.15) 

where 1/A, is the average rate of exponential growth of the solution. There is increased interest in \t | 2 because it is related 
to the conductance or resistance of the system. Landauer considered a disordered segment of length L connected to 
2 semi-infinite leads. The leads are considered to be connected to reservoirs at ±°° and all relaxation of energy takes 
place in these reservoirs. Then we can write equations for the right and left going contributions to the current on both 
sides of the system 

ijw. = |*| 2 »£f + M 2 »*- 

, , (8-16) 

i R - = \r\ z i L+ + \t\ z i R . 

where the letters L and R refer to the left and right side of the disordered region respectively and the signs + and - to 
the direction of the current. 

After relaxation to the ground state the change in potential energy is proportional to the total occupancy of the 
states which is proportional to the scalar sum of the 2 currents. The number of states associated with a change of 
potential energy ebV is 

1 dk 

N s = eZV- YE (8.17) 
whereas the absolute current carried by each state is calculated using the group velocity 

1 dE 

l * = %Tk (8 ' 18) 

Thus the total current associated with 8V is the product of these 

e 1 

i = NJ s = 2—8V (8.19) 
h 

where the factor 2 is for spin. 

We consider only right flowing total current (i.e. z' r - = 0) and combine (8.16) with (8.19) to obtain 

e 2 |f| 2 

G=2 (8.20) 

h 1 - |f| 2 
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This formula has certain attractions. It does all the right things in the limits. At zero disorder |f| 2 >— >1 and G is infinite. 
However the formula is not uncontroversial. Other authors (e.g. J.B.Pendry) prefer the form 

G = 2— \t\ 2 (8.21) 
h 

which is obtained by ignoring the contribution of the reflected current to the voltage difference. The differences 
between these formulae really lie in the description of the leads rather than the disordered region itself. For example 
consider the case of zero disorder. (8.21) would allow the potential to drop solely over the non-existent disordered 
region: clearly absurd. On the other hand in the case of a narrow channel between wide leads there is no difficulty. 

The root cause of this controversy lies in a fundamental contradiction. Conductivity implies transport between 
2 points at different electrochemical potentials. Hence transport and dissipation are inextricably linked. On the other 
hand there clearly is a relationship between transmission coefficient and conductance. Where does the dissipation 
come in? 

In addition to the above uncertainties in defining properly the DC resistance of a finite one-dimensional system, 
it is found that the probability distribution of |f| 2 possesses long tails which are responsible for a peculiar behaviour of 
various averages. For example 

(ln|f| 2 ) = -2L/X (8.22a) 
(|ff 2 ) = \ [1 + exp(4ZA)] (8.22b) 
(|f|" 4 ) = i + \ exp(4ZA) + \ exp(12ZA) (8.22c) 

(|f| 2 ") D n (kl2L) m exp(-L/2A.) (8.22J) 

where the constant D„ depends only on n. The standard deviation of In \t | 2 is proportional to (LIX) 1/2 for large LIX, with 
a distribution that appears to be Gaussian. 

The localisation length can be obtained from the Greens function as follows 

AT 1 =-— !— -(\n\G(l,m)\) \l-m\^oo (8.23) 
\l — m\ 

Using the Renormalised Perturbation Expansion (see earlier) we have 

G(l, m) = G{1, l)VG(l +1,1+1 [l])V . . . G(m, m[m - 1]) (8.24) 

Note that in 1-d G(m,m[m - 1]) is independent of everything to the left of site m. It can be shown that products 
like (8.24) have the property that, if we define the submatrix G' to contain all the matrix elements G(z'j') with ;' and j 
between I and m, then 

G(/,m) = det{G'} (8.25) 



Hence (8.23) becomes 



lim /y^ln I £' — £',• I 

|Z-m|-*» \l-m\ \ ^ 1 71 



/ E,5(F- g ,U , . (8 . 26) 



= lim / dE' -)ln\E-E' 

|/-m|->oo J_^ \ \l — m 



= J (p(E')) In \E-E'\ dE' 
This result was first obtained by Thouless. 



b) Scaling Theory 

The major breakthrough came in 1979 with the scaling theory of Abrahams, Anderson, Licciardello and Ramakrishnan 
(the gang of four). Many of the ideas go back to Thouless and to Wegner. 
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i) What's so special about the Conductance? 

Consider a block of disordered material of dimensions L d . The average difference in energy between neighbouring 
eigenstates is AE = L d lp where p is the density of states. Another important energy scale is the effective uncertainty 
in the energy of each state, 8E. When the disordered block is embedded in a solid we can define 8E in terms of the 
average time a particle spends in the block. Thus 

SE = h/x (8.27) 

When L is larger than the mean free path the motion is effectively diffusive. Hence x is related to the size of the block 
via the diffusion constant, D. 

L 2 = Dx (8.28) 
D however is related to the conductivity via the Einstein relation 

a = e 2 Dp (8.29) 



Hence we can write 

8E = %op- l L- 2 = %aAEL d ~ 2 = gAE. (8.30) 

e l e l 

The physical picture is that if hE <C AE electrons from one block will not penetrate into the next block when several 
blocks are joined together, whereas if 8E > AE there will be no barrier to transport. However (8.30) shows that 8E/AE 
is essentially the conductance and this is the only quantity which controls the behaviour as blocks are joined together. 
This result can be summarised in the form 



g(2L)=f[g(L)] (8.31) 

or in differential form 

4 ' =pfe) or ^M = pOng) (8.32) 



dlnL ro dlnL 

This is the famous p function. Two important points should be borne in mind: 

a) (8.32) is simply a special case of a more general multi-dimensional relationship. 

b) A more general form would be to consider the whole distribution of values of g rather than an unspecified mean. 
So that we could write 

d P(lns) = p[P(lng),lng] (8.33) 



dlnL 

We shall return to this case later. 



ii) General Scaling Properties. 

(8.32) looks like a non-linear equation of motion such as is studied in the theory of chaos, except that ft— > InL. Most 
of the mathematical results from chaos can be applied to our problem. In particular if we consider a multi- or infinite- 
dimensional version of (8.32) it will in general be possible to identify various topological features of the phase space, 
such as 

a) Attractive and repulsive lines, surfaces, etc., where the flow is always towards (away from) the feature. 

b) Fixed points. These are of 3 types. 

1) Stable fixed points or sinks. Everything flows towards the point. 

2) Unstable fixed points or sources. Everything flows away from the fixed point. 

3) Saddle points. Flow is towards the point in some directions, away in others. 

It should be noted that it is always possible to write a general equation like (8.32) but it may have to be infinite 
dimensional. The one-parameter scaling theory contends that there exists an attractive line which may be parameterised 
by the conductance g. Thus we must be careful to make sure that the system is large enough that it is sufficiently close 
to that attractive line for the assumption that g is the only relevant variable to be valid. The other variables which 
control the deviation from the attractive line are often termed irrelevant. In reality this does not necessarily imply that 
they may be ignored. 
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iii) Fixed Points 

Close to a fixed point an equation like (8.32) may be linearised to give 



d? =B(g-g*) (8.34) 



dlnL 

where B is the matrix of derivatives of p. (8.34) has solutions of the form 

g i = g*+A i L^ (8.35) 

where the p', are the eigenvalues of B. The constant of integration A, contains all the information about disorder, 
energy etc.. Consider for example one of the solutions (8.35). For any finite L we expect g to be a smooth function of 
disorder, 8e, or energy, E. In this case we may Taylor expand A around the fixed point. Hence 

g = g*+A'(E-E*)LV (8.36) 

where E* is the fixed point energy, i.e. the energy at which g is independent of L. It is useful to write everything in 
terms of dimensionless variables. Thus 

g = g*+A"(Ll^f (8.37) 
Comparing (8.36,37,38) we obtain an expression for the length scale £ 

^~ \E-E*\~ V V (8.38) 
from which we can identify the critical exponent v, associated with the correlation length, as 

Vj = -1/p'j (8.38) 

iv) Crossed Diagrams 

From the totally crossed diagrams (see beyond the CPA above) we can obtain an expression for the expansion of p(lng) 
in terms of llg. Suppose the leading terms (in 3-d) are 

££=.-(£)'- »(W-.n S -) 

which leads directly to the result that 

v = lln (8.40) 

Thus if the leading term in the expansion is llg then V = 1, whereas if the leading term is llg 1 then v = \. Note 
however that this is a very rash approximation. Recent results (Wegner) suggest that the next term in this expansion is 
at least as large as this. Nevertheless for ordinary disorder we have V = 1, but for random magnetic field or spin-orbit 
coupling v - 1 



2" 



v) Conductivity 

For large g we expect the scaling behaviour to be classical, 



g = ^oL d ~ 2 (8.41) 



which gives directly the limiting behaviour p = d - 2. However, comparing (8.41) with (8.37) we obtain o ~ and 
hence the critical exponent s for the conductivity is related to V by s = (d - 2)v. This result, originally derived by 
Wegner, was often misinterpreted to predict a minimum metallic conductivity in 2-d (s = 0). However it presupposes 
the existence of a fixed point, which is not necessarily valid in 2-d. 
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vi) Whole Distribution 

Returning to the general (8.33) for a one-parameter scaling relationship for the whole distribution of g. If there exists 
a fixed point here it is represented by a scale invariant form of the distribution. In this case the linearisation takes the 
form 

^-= B [P(g)-P t (g)] (8.42) 

which has the solution 

P(g) = P*(g) + A(g)L B (8.43) 

where A(g) depends on energy, disorder etc. and A(g) dg = 0. Note, in particular, that the value of the critical 
exponent is independent of which average is considered, (g), exp (\ng) or 1/ (II g), as long as the corresponding integral 
over P(g) is well defined. It is by no means obvious that this last condition is fulfilled. 

Another side effect of the existence of a fixed point distribution is that all the moments of P*(g) are also fixed. It 
may well appear that these moments are independent of the size of the system. There is considerable evidence for the 
existence of a universal value for the standard deviation of the conductance. 

e 2 

= 1 or 8G = — (8.44) 
n 

These are often called Universal Conductance Fluctuations. The effect is observed at very low temperatures by varying 
Fermi level (gate voltage) or magnetic fields. 

c) Field Theory Formulation 

This section can only be the briefest of outlines of the field theoretical formulation of the problem. 
Consider the integral form 
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D0 1 D0 exp {i «>t [E + iy- e] 0) } = ^ +iy _^ (8.45) 



Note that this is really a 2-dimensional integral, because we have integrated over the and T separately. The more 
general form of this is then 

/ + °° (\k) n 
_ UD^D^j expji (J>T [E + iy-Hy] 0,)} = - ^ + [ y _ H] (8.46) 

From this form it is possible to calculate the Green's function using 

J + fp0fD0 ; 0t0„exp{...} = G m „^^ (8.47) 

or even the 2-particle Green's function by including 4 (j)'s in (8.47). However the presence of the determinant (8.47) 
makes it difficult to evaluate exactly. There are 2 ways to solve this problem: 

i) The Replica Trick. Consider n replicas of the system. Then the total system will include (8.46,47) raised to the 
power n. If we then take the limit n — > we obtain det [.. .]° = 1 and are left with the desired Green's function. 

ii) Supersymmetry. When the integral in (8.46,47) is evaluated using Grassmann or anticommuting variables, the 
result is the same except that the determinant appears in the numerator. By integrating over 2 normal (commuting) 
variables and 2 Grassmann variables the determinants cancel and again we are left with the Green's function as 
required. 

In the case where the 2-particle Green's function is required (e.g. in order to calculate the conductivity), there is 
a further complication. The integral in (8.45,46,47) only converges when the real part of the exponent is negative. 
In (8.45) a factor i is introduced which combines with the imaginary part of the energy iy to give -y0 2 . For the 
conductivity, however, we require terms like G + G~. To enable us to calculate such quantities the integrand must take 
the form 

i U\. [E + iy-H^hj- 2/ [E - iy- H y ] 2i ) 

f+ ♦ 1 f ■ t 1 (8 ' 48) 

= i {nMj - [E-Hij] + iy jtf^y + 02,42/ j 
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The value of the exponent in the integrand in (8.48) is made up of a term which is invariant under transformations 
which preserve the value of ((j)};(|)i/ - §l$2j), and a small perturbation, proportional to y which breaks this symmetry. 
The group of allowed transformations is the pseudo-unitary group, U(n,n). In the usual case where H is real and 
symmetric the real and imaginary parts of (j) are independent of one another and the group is the pseudo-orthogonal, 
0(2n, 2n). When the only source of scattering is the spin-orbit coupling there is an extra symmetry due to the spin 
exchange under time-reversal symmetry. The resulting group is pseudo-symplectic, S(n, n). 

These 3 cases form the Universality Classes for the localisation problem. In 3-d the orthogonal case, ordinary 
scattering, gives a critical exponent of 1, whereas the Unitary case gives v = s = i. There are no reliable predictions 
for the Symplectic case (the series doesn't converge). Experimentally this is still a matter of some controversy. In 
heavily doped semiconductors, for example, the exponent appears to be unity or \ for compensated and uncompensated 
systems respectively. 

In 2-d the 3 cases differ in the way the P-function approaches zero, giving p = — alg, p = -alg 1 and P = +alg 
respectively. Note in particular the positive contribution from the symplectic (spin-orbit) case. This effect is sometimes 
termed, weak antilocalisation. Note that this implies the existence of extended states. It is not completely clear whether 
in fact there are any localised states in this case. 

Note however that in real systems there will be a variety of scattering mechanisms, such that no single universality 
class need apply. 

d) Numerical Approach. 

A number of numerical approaches to the localisation problem have been developed. In the 1970's there was a 
remarkable degree of agreement on the details of the transition in 2-d. In retrospect the agreement was spurious, being 
due to the similar storage capacities of the various computers used. 

More recently finite sized scaling methods have come to be accepted more generally. Using a numerical variant 
of the transfer matrix technique 

'y n +i\ _ ( E-H n -A / y„ 
y n J \ I J \y n -! 



-u 

m=l 



E-H m -A y x 



(8.49) 



= T n 



we can calculate the Lyapunov exponents, essentially the logarithm of the eigenvalues of T. These are in pairs which 
are reciprocals of each other due to the symplectic symmetry of T. The smallest exponent may be identified with the 
localisation length X M of a strip or bar of width M(xM) for which the calculation is being carried out. 

It is found (MacKinnon & Kramer) that A = Xm/M obeys a similar scaling relationship to the conductance. Thus 
many of the results discussed above can be transferred to A. In particular, by fitting 

A = A* +A(W— W*)M a (8.50) 

to results for A as a function of strip width M and disorder W, it is possible to calculate the critical exponent. Most 
work has been done on the Rectangular distribution which seems to give s = v » 1.5 which is not compatible with 
the analytical approximations or with the experiment. However more recent work using other distributions is in better 
agreement with analytical work s = V = 1 . 

Clearly however the problem is not yet solved. 
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